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J  Abstract— This  paper  giiev  several  basic  results  on  dynamic  nonlinear 
__  networks  from  a  geometric  point  of  view.  One  of  the  main  adtantages  of  a 

\  geometric  approach  is  that  it  is  coordinate-free,  i.e.,  results  obtained  by  a 

c  geometric  method  do  not  depend  on  the  particular  choices  of  a  free,  a  loop 
pH  matrix,  state  variables,  etc.  Therefore,  the  method  is  suitable  for  studying 
intrinsic  properties  of  networks. 

Cgf*  It  is  shown  that  transversality  of  resistor  constitutive  relations  and 
Kirchhoff  space  is  a  sufficient  condition  for  the  configuration  space  to  be  a 
submanifold.  Main  result  of  the  paper  states  that  a  network  is  locally 
,  f  solvable,  i.e..  the  dynamics  of  a  network  is  well  defined  in  the  sense  of 
Definition  3.  if  and  only  if.  capacitor  charges  and  inductor  fluxes  serve  as  a 
local  coordinate  system  for  the  configuration  space.  In  other  words,  if  all  the 
variables  in  a  network  are  determined  as  functions  of  capacitor  charges  and 
inductor  fluxes,  at  least  locally,  then  the  dynamics  is  well  defined.  Con¬ 
versely.  if  the  dynamics  is  well  defined,  then  all  the  variables  in  a  network 
are  determined  as  functions  of  capacitor  charges  and  inductor  fluxes. 
Because  of  its  coordinate-free  property,  the  main  result  also  says  that  if  the 
dy  namics  is  well  defined  in  terms  of  some  coordinate  system,  then  it  must 
be  well  defined  in  terms  of  capacitor  charges  and  inductor  fluxes.  Con¬ 
versely.  if  the  dynamics  is  not  well-defined  in  terms  of  capacitor  charges 
and  inductor  fluxes,  then  there  is  no  choice  of  variables  in  terms  of  which 
the  dynamics  is  well  defined  in  the  sense  of  Definition  3.  This  shows  that 
capacitor  charges  and  inductor  fluxes  are  the  fundamental  quantities  in 
describing  the  dy  namics  of  networks.  Perturbation  results  are  given  which 
guarantee  transversality  and  local  solvability.  Finally,  several  other  per¬ 
turbation  results  are  given  which  guarantee  eventual  strict  passivity  of 
dynamic  nonlinear  networks.  They  explain  why  the  voltage  and  current 
waveforms  of  almost  even  network  of  practical  importance  are  eventually 


uniformly  bounded^ 

/ 1*  ,  r  »  . 

,  f  +  /fy  ( 


i*  ntrodi  ction  means  that  if  all  the  variables  in  a  network  are  determined 

THIS  PAPER  gives  several  basic  results  on  dynamic  as  functions  of  capacitor  charges  and  inductor  fluxes,  at 
nonlinear  networks  from  a  geometric  point  of  view,  least  locally  at  each  point,  then  the  network  is  locally 
One  of  the  main  advantages  of  a  geometric  approach  is  solvable,  i.e..  the  dynamics  is  well  defined.  Conversely,  if  a 
that  it  is  coordinate-free,  i.e.,  the  results  obtained  by  a  network  is  locally  solvable,  then  all  the  variables  in  the 
geometric  method  do  not  depend  on  the  particular  choices  network  are  necessarily  determined  as  functions  of  capaci- 
of  a  tree,  a  loop  matrix,  state  variables,  etc.  Also,  the  tor  charges  and  inductor  fluxes.  Because  of  its  coordinate- 


geometric  approach  allows  us  to  choose  a  convenient  coor¬ 
dinate  system  and  use  it  to  derive  general  conclusions 
which  hold  with  respect  to  any  other  coordinate  system. 
Therefore,  this  approach  is  suitable  for  studying  intrinsic 
properties  of  networks  and  it  enables  us  to  resolve  and 
clarify  a  number  of  subtle  paradoxes  and  perplexing  ques¬ 
tions  which  lie  at  the  very  foundation  of  nonlinear  circuit 
theory.  In  particular,  several  basic  questions  involving  the 
formulation  of  state  equations  for  nonlinear  networks  are 
hereby  resolved  in  a  rigorous  manner. 

In  Section  II  we  will  describe  nonlinear  networks  in  a 
coordinate-free  manner.  In  Section  III  we  will  discuss 
transversality  of  the  resistor  constitutive  relations  and  the 
Kirchhoff  space.  Transversality  is  important  in  that  it 
guarantees  the  configuration  space  to  be  a  submanifold. 
We  will  give  two  perturbation  results  which  guarantee 
transversality.  One  involves  element  perturbation,  i.e..  per¬ 
turbing  the  existing  resistor  constitutive  relations.  The  other 
involves  network  perturbations,  i.e..  augmenting  the  net¬ 
work  with  capacitors  and  inductors.  In  Section  IV  we  will 
discuss  local  solvability  which  is  a  condition  for  the  dy¬ 
namics  to  be  well  defined.  Main  result  ( Theorem  1)  says 
that  a  network  is  locally  solvable,  i.e..  the  dynamics  is  well 
defined  in  the  sense  of  Definition  .3.  if  and  only  if.  capacitor 
charges  and  inductor  fluxes  sene  as  a  local  coordinate 
system  for  the  configuration  space  of  the  network.  This 
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free  property,  the  main  result  also  says  that  if  the  network 
is  locally  solvable  in  terms  of  some  coordinate  system,  then 
it  must  be  locally  solvable  in  terms  of  capacitor  charges 
and  inductor  fluxes.  Conversely,  if  the  network  is  not 
locally  solvable  in  terms  of  capacitor  charges  and  inductor 
fluxes,  then  there  is  no  choice  of  variables  in  the  network 
in  terms  of  which  it  is  locally  solvable.  This  says  that 
capacitor  charges  and  inductor  fluxes  are  the  fundamental 
quantities  in  describing  the  dynamics  of  networks.  One  of 
the  interesting  implications  ( Corollary  5)  of  this  result  is 
that  if  a  network  is  locally  solvable,  then  the  capacitors  are 
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necessarily  locally  charge  controlled  and  inductors  are  nec¬ 
essarily  locally  flux  controlled.  This  seems  to  explain  why 
almost  every  capacitor  (resp.  inductor)  of  practical  impor¬ 
tance  is  locally  charge  (resp.  flux)  controlled.  After  proving 
the  main  result,  we  will  give  a  network  perturbation  tech¬ 
nique  which  guarantees  local  solvability.  In  Section  V  we 
will  give  several  perturbation  results  which  guarantee  even¬ 
tual  strict  passivity  of  dynamic  nonlinear  networks.  Theo¬ 
rem  2.  another  important  result  of  the  paper,  explains  why 
the  voltage  and  current  wave  forms  of  almost  every  net¬ 
work  of  practical  importance  are  eventually  uniformly 
bounded. 

General  Remark:  For  simplicity,  we  will  usually  delete 
the  superscript  T  denoting  the  “transpose”  of  a  vector  or 
matrix  whenever  no  confusion  arises. 

II.  Coordinate-Free  Description  of  Nonlinear 
Networks 

Throughout  the  paper,  we  need  to  use  the  fact  that 
transversality,  local  solvability  and  eventual  passivity  arc 
coordinate-free  properties,  i.e..  they  are  independent  of  the 
choices  of  a  tree,  a  loop  matrix,  a  cut  set  matrix,  state 
variables,  etc.  Here  we  will  explain  how  nonlinear  networks 
are  described  in  a  coordinate- free  manner. 

Consider  a  nonlinear  network  Ai  containing  nK  resistors. 
nc  capacitors  and  n,  inductors.  Let  r  ~H'H.  v( .  ty )  and 
i=(i„.  i(  ,  /,  )  be  the  branch  voltages  and  branch  currents, 
respectively,  and  let  q  and  <j>  be  the  capacitor  charges  and 
inductor  fluxes,  respectively,  where  R.  C.  and  L  denote 
resistors,  capacitors,  and  inductors,  respectively.  The  fol¬ 
lowing  are  the  standing  assumptions  of  this  paper: 

(a)  The  linear  graph  0  which  defines  the  topology  of  At 
is  connected. 

(b)  At  is  time  invariant. 

(c)  The  resistor  constitutive  relations  are  characterized 
by 

(r,.i,)6A,CRJ**  (I) 


Connectedness  is  necessary  for  tree  to  exist. 

2)  Most  of  the  results  of  this  paper  can  be  easily  gener¬ 
alized  to  include  the  lime-varying  case  under  appropriate 
conditions.  We  make  this  assumption  simple  to  avoid 
introducing  complicated  notations. 

3)  Under  assumption  (c)  resistors  can  be  coupled  to  each 
other  and  they  need  not  be  voltage  or  current  controlled. 
This  includes  virtually  all  modes  of  representation,  includ¬ 
ing  the  hybrid  and  transmission  representations.  In  partic¬ 
ular.  a  broad  class  of  nonlinear  dependent  sources  arc 
allowed  in  this  formulation.  We  regard  independent  sources 
as  uncoupled  resistors.  All  multi-terminal  elements  are 
represented  as  coupled  two-terminal  elements. 

4)  Under  the  present  formulation,  capacitors  need  not  be 
voltage  or  charge  controlled.  Similarly,  inductors  need  not 
be  current  or  flux  controlled.  Notice  that  capacitors  can  be 
coupled  to  each  other  and  inductors  can  be  coupled  to  each 
other. 

5)  Coupling  among  elements  of  different  kinds  are  not 
allowed.  For  example,  dependent  sources  controlled  b\ 
variables  of  reactive  elements  are  not  allowed. 

6)  We  need  C 2  property  of  A*.  A( .  and  A,  rather  than 
C'1  because  we  would  like  to  define  C1  vector  fields  on  the 
configuration  space.  (See  Section  IV.) 

7)  Assumption  (e)  was  introduced  only  for  simplicity. 
This  involves  no  loss  of  generality  in  view  of  the  results  of 
Chua  and  Green  [1]  and  Sangiovani-Vincentelli  and  Wang 
(21- 

Now  let  h nK  +  nc+nl  and  let 
A-  {(v.i.q.<l>)\(t'R.iR)e\H.(i\  ,q)E\i  .(r,  A,  }. 

(6) 

Then  it  follows  front  (c)  and  (d)  that  A  is  a  (/>  *  n,  t  «,  )- 
dimensional  C:  submanifold.  Let 

K-  {( r. /.</,</>)!( l'-  *)  satisfies  Kirchhoff  Laws] .  (7) 


where  \H  is  an  /t„ -dimensional  C2  submanifold, 

(d)  Capacitors  are  characterized  by 

(vc.q)e.\cCRln‘ 


and 


<iq 

dt 


c 


(2) 

(3) 


where  A,  is  an  n(  -dimensional  C2  submanifold.  Inductors 
are  characterized  by 


and 


A,  CRJ"' 

d$  _ 

dt  V' 


(4) 

(3) 


where  A,  is  an  n ,  -dimensional  C:  submanifold. 

(e)  There  are  no  capacitor-only  loops  and  no  inductor- 
only  cut  sets. 

Remarks:  1)  There  is  no  loss  of  generality  in  assuming 
(a)  since  disconnected  subgraphs  can  be  hinged  together. 


It  is  well  known  that  K  is  a  (/>  (»,  f  nt  (-dimensional 
linear  subspace.  This  space  is  called  the  Kirchhoff  space  and 
is  independent  of  the  particular  choices  of  a  tree,  a  htop 
matrix,  a  cut  set  matrix  etc.  Since  ( t\  i.q.  $ )  must  satisfy 
the  constitutive  relations  and  the  Kirchhoff  laws  simulta¬ 
neously.  the  operating  points  are  restricted  to  within  the 
follow  ing  subset  . 

2=  A  O  A .  (8) 

The  set  -  is  called  the  configuration  space  of  Ai  since  this  is 
where  the  dynamics  takes  place. 

III.  Transyirsai  it  y 

Since  the  dynamics  takes  place  on  the  configuration 
space  2!.  the  object  -  should  be  well  behaved  at  least  to  the 
extent  that  we  can  write  down  differential  equations  on  it 
For  that  purpose,  it  suffices  to  require  2  to  be  a  differen¬ 
tiable  submanifold.  A  little  problem  is  that  even  if  A  and  K 
are  perfectly  well-defined  differentiable  submanifolds,  their 
intersection  -  may  or  may  not  be  another  submanifold.  A 
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sufficient  condition  for  2)  to  he  a  submanifold  is  the  Proof:  It  follows  from  a  similar  argument  to  the  proof 

transversal! ty  [3]  of  A  and  A.  which  is  abbreviated  by  of  Proposition  I  of  (3|  that  AffiA  if  and  only  if  for  each 
A  A .  It  is  shown  in  [4]  that  if  AffiA.  then  2)  is  an  arG^. 

(n<  +  nt  )-dimensional  submanifold.  This  is  true  for  any  C  «'  i  q  <t> 

sub  ranifolds,  r>  1.  We  first  give  a  formula  for  checking  \  B 

transversality  of  A  and  A.  .. 

Since  A  „  is  a  C2  submanifold  of  dimension  nR,  for  each  ^  v 

point  (r„  ti,  iR  )G A*.  there  is  a  neighborhood  UR  CR:"*  of  fn 

this  point  and  there  is  a  C2  function  fR\  UR  — ♦R"*  such  :D.f(  D  f 

that  r‘  * 

:  d„  f,  D+f, 


A KnuK=fR  '(0) 


rank (/>/*)„.  ,  =  nR.  for  all  (t,.i,)6A,n  UR  where  BJS  as  *n  <15)  7and  ,s  the  fundamental  cut  set 
R  "*•'**  R  R  matrix.  Since  Q  =  \-Bl,  :  1),  one  can  show,  by  elemen- 

(10)  tary  operations,  that  (18)  holds  if  and  only  if  (16)  holds.  □ 

where  (Dfn\r  ,  ,  is  the  derivative  of  fR  at  <r R.iK).  Simi-  Remark:  "T  need  no1  be  a  ProPer  lree-  ()ne  simP'>  has 

larlv.  for  each  point  (t;  .?„)6A(  (resp.  {i,  ,  <j>0  )£.\( ).  !wo  more  nonzcro  subma.rices  in  (17). 

there  is  a  neighborhood  t/"cR2"*  (resp.  U,  CR;"' )  of  this  lf  A<  A,  admit  special  forms,  then  we  can  give 


point  and  there  is  a  C2  function  f(  :  U(  — »R',<  (resp .  f,  : 
V,  — ►  R"' )  such  that 

A,  D  i/\-  =  fc  ‘(0)(resp.  A ,  D  V,  -f,  '(0))  (11) 

and 


more  explicit  formulas. 

Definition  l:  Capacitor  constitutive  relations  A,  is  said 
to  be  locally  voltage  (resp.  charge)  controlled  if 

rank (0, =nr,  for  all  ( t\.q  )E  A(  nl' 

(19) 


rank(D/.)  =  n,.,  for  all  ( a  .q )  G  A(  D  U.-  (12)  ,  .  .  „ ,  .  ,  , 

*"*•*'  '  (resp.rank  (D,.fc)(=nc.  for  all  (c<  .q)e A<  n  L[  ) 

(resp.  rank  (Df,\  ,+,  =  n, .  for  all  (i,..^)€  A,  O  U, ).  ....  -  '.  ..  , 

where  U(  is  as  in  (11)  and  (12).  Similarly,  inductor  con¬ 
it  follows  from  (e)  that  there  is  a  proper  tree  ‘.'T.  Let  C  be  stitutive  relations  A,  are  said  to  be  locally  current  (resp. 
its  associated  cotree  and  let  c  and  i  be  partitioned  accord-  flux)  controlled  if 

ingly:  ,  ,  4  rank  (/>,/,),.  *=«,.  for  all  (/,.*)€  A ,  f)U, 

r=(tv  :  e.;)=  .t>,  :  vR  .  er  (13)  ,  ,  ,  . 

'  ( resp .  rank ( D,;  f, ), (/  ^ ,  =  n , .  for  all  («,  ,<J*)e  A,  n  V, ) 

/  (/,.  .  )  (,/v,7  •  'n,  •'()■  ('4)  where  U,  is  as  in  ( 1 1 )  and  ( 12). 

Let  B  be  the  fundamental  loop  matrix  associated  with  VL  Let  tt'R:  R;,M  -»IR  2"s  be  the  projection  map  defined 


»  I1  *1 

Set  x-  ( e.  i.q.  <f>  )- 

Proposition  I:  Affi A  if  and  only  if 

rank^(ar)  =  b,  for  all  *G  2 

where 
9(x)  = 


(15)  «R(c.i\q.<j>)-(vR.iR).  (20) 

Let  i:  --*ft2h'"‘  be  the  inclusion  map  defined  by 

(v.i.q.Q)  (21) 

(16)  and  set 

•nR--nR°i.  (22) 

We  next  decompose  B  of  (15)  according  to  (13)  and  (14): 

B i  -  Brh  Br(  .  (23) 

®/.«  ®/i 

Corollary  I:  l  et  A(  (resp.  A,  )  be  locally  voltage  (resp. 
current)  controlled.  Then  AffiA  if  and  only  if 

'2*  TanY.9R(vR.iR)  =  nR.  for  all  ( t>.  »„  )£«/,( Z)  (24) 


Here  Drii  fR  denotes  partial  derivative  of  fR  with  respect  to  Proof:  It  follows  from  (17)  that  if  A,  (resp.  A,)  is 
r and  denotes  a  zero  submalrix.  Other  symbols  have  locally  voltage  (resp.  current)  controlled,  then  (16)  holds  if 
similar  meanings.  and  only  if 


DrJR  ~{DCtfR)B,  ; 

oitfK  +  (o„/J*r: 

;  Drf( 

»,/< 

•  ;•  Di,  ft  : 

(17 

SilUVi  >»„)  = 
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rank  [  DrJR 

-(ft*/j.)ft 

:  fte/«+ (ft, /*)«'],,,,, 

(26) 

Since 

ft,  f* 

[  D*J*  :  • 

]•  ft*/*  -[ft,/*:  ■] 

ft,/«  = 

;  ] 

■  i>,Jh  =  [  d,Jr  :  •] 

substituting  these  and  (23)  into  (26),  we  have  (25).  Since 
( i\  i )  £  -.  the  vector  ( vR.  iR )  must  belong  to  ir^  2 ).  □ 

Corollary  2:  Let  A,  (resp.  A, )  be  locally  charge  (resp. 
flux)  controlled.  Then  AffiA  if  and  only  if 

rank^f2(x) =  nR,  for  all  (27) 

where 


1  ig  I  A  nonlinear  circuit  »nh  A  <T.  A  (a)  I  lie  circuit  diagram  (h) 
C  apacitor  constitutive  relation  is  locall\  \oltugc  controlled  <c)  C'apaci 
lor  constitute c  relation  is  local K  charge  controlled 


(resp.  current)  controlled.  Then 

[(«:  {OF)y:)  (a,  ( />A  )y  l)BRR  : 

(«,  (DF)yt)BRl  :  O,  (DHS,) 

+  (ft  (DF)6:)BfR  :  </J,  (Of  )«,)*/ i(ii 

(31) 


~{DriiJR)BRR  :  DiK  fR  t{D,JR)BZK  ;  (».  „/«)*/«  (».  ,  /,  )  :  (».„,/* )*/„(!>,,/,  )  '/>*/,]  t-  (28) 


Proof:  If  A(  (resp.  A;  )  is  locally  charge  (resp.  flux) 
controlled,  then  ( £>,  /,  )  (resp.  ( /),  /,  ))  is  nonsingular. 
Therefore,  bv  elementary  operations,  one  can  show  that 
(18)  holds  if  and  only  if  (27)  holds.  U 

Example  I :  Consider  the  circuit  of  Fig.  1(a).  where  A, 
is  described  by  q=gif 1\  )  as  in  Fig.  1(b).  Then  V?  -  {C ' }  is  a 
proper  tree.  BRC  -  1  and 

•'«( '«''«)  =  [  1),Jk  ;  t>, 

which  has  rank  I  because  of  ( 10).  Therefore.  A  A  for  any 
A  R  as  long  as  it  is  a  (  2  submanifold.  Suppose,  now.  that 
i\  -gc(q)  as  in  Fig.  1(c).  Since  g,  is  not  injective.  A,  is 
not  locally  voltage  controlled  and  Corollary  I  does  not 
apply.  In  order  to  apply  Corollary  2,  we  compute 

'f:u>  [/>„/« ;  (",«/*)( />«< )], 

which  may  or  may  not  have  rank  1  depending  on  fR  and 
X,  .  If  iR  --/„(  rR ).  however,  the  above  matrix  always  has 
rank  I  and  A  ft  A". 

Next  suppose  that  A  R  admits  a  generalized  port  coordi¬ 
nate  (3).  i.e..  \R  is  represented  by 


where  a.  P .  y,  and  8  are  nR  n R  matrices  and  F  .'  R""  — » 
R""  is  a  (  ’  function.  Recall  the  partition  t-R  -( r„f  :  vR^ ), 
iH  {iR  iR , )  and  partition  a.  p,  y,  and  8  accordingly: 

<*-[<»,  :  or,].  P  [p,  ft] 

Y  [yi  ■  Y:]-  «■  [«|  :  ft]  (30) 

Also  recall  that  A  R  is  said  to  be  globally  collage  controlled 
(3|if{  iR,  i)  vR  and  globally  current  controlled  if  £  vR.  tj 

*V 

Corollary  3.  Let  A  R  admit  a  generalized  port  coordinate 
representation  and  let  A,  (resp.  A,)  be  locally  voltage 


In  particular,  if  \R  is  globally  voltage  controlled,  then 


(32) 

and  if  A  R  is  globally  current  controlled,  then 

BRR]  ««,  I 

:  ;  </>n  f  : 

1  ’  *  K 

inn 

k«ll,  W) 

Wc  can  also  give  similar  formulas  for  9'  when  A,  (resp. 
A,  )  is  locally  charge  (resp.  flux)  controlled.  We  omit  them, 
however. 

Now  suppose  that  A„  is  globally  parametnzable  |3|.  i.e.. 
A  R  is  globally  diffcomorphic  to  R"»  and  write 

(  1'k(Ph  )•  *«(  Pm  ))  ^«(Pk)'  Pm  1  -R  (33) 

where  A  R  —> ►  R"*  is  a  global  coordinate  system. 

Definition  2:  A,  (resp.  A,  )  is  said  to  be  globally  para- 
metrizuble  if  (  r,  ,  q ) t£  A,  ( resp.  ( i, ,  ft) e  A ,  )  is  represented 
by 

(*'<  (P<  )-V(ft  ))  *<  '(ft)-  P<  t:R"'  (35) 

( resp .  (i,(p,  ))  '(p,  ).p,  CR"> ) 

where  if(  :  \(  — >R'''  (resp.  \p, :  \,  — ►  R ” ' )  is  a  global 
coordinate  system. 

If  \  R.  A,  and  A,  are  globally  parametrized,  then  A  of 
(6)  is  parametrized  by  p  (p„.p,  .  p,  )  and  (i,  .  r,  ).  We 
write  this  as  x  .r(p.  i(  .  r,  )t  A. 
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Proposition  2:  Let  A „.  A,  ,  and  A,  he  globally  parame¬ 
trized.  Let  "5  be  a  proper  tree  and  let  >:  be  its  associated 
cotree.  Then  A ffiA  if  and  only  if  for  each  p€ ft with 
x(p.  /( .v,  )e2: 

rankf*(p)  =  «w  (36) 


f*(p) 


B*r(  Dt'ft, 

Brr(  D‘k(  i 


Proof:  First,  observe  that  for  any  .r€.\ 

Dv 

Di 


Recall  that  the  Kirehhoff  space  A  is  parametrized  by 


v  =  Q'i 

This  implies  that 


qew"* .  <>enr'. 


T.K  =  Im 


By  definition  [3).  A  ft  A  if  and  only  if 


Dq 

.r,i 

which  is  equivalent  to 


J  J  ( n  i,  i 

jatrix  has  the  following  form: 


By  elementary  operations,  one  can  show  that  (41)  holds  if 
and  only  if  (36)  holds.  L 

Remark:  Proposition  2  holds  even  when  A  w.  A,  .  and  A , 
are  locally  parametrized  at  each  point,  and  it  includes  [5]  as 
a  special  case. 

Suppose  now  that  AiJfA'.  Then  it  would  be  helpful  if  one- 
can  perturb  ^  in  an  appropriate  way  such  that  the  result¬ 
ing  network  satisfies  transversality.  In  the  following  we 
give  two  perturbation  results.  The  first  method  involves 
element  perturbation  and  consists  of  perturbing  the  exist¬ 
ing  resistor  constitutive  relations  A#.  The  second  method 
involves  network  perturbation  and  consists  of  augmenting 
-At  bv  adding  arbitrarily  small  linear  inductors  and  arbi¬ 
trarily  large  linear  capacitors  by  pliers-type  entry,  and  bv 
adding  arbitrarily  large  linear  inductors  and  arbitrarily 
small  linear  capacitors  bv  soldering-iron  entry.  Therefore, 
in  the  limit  we  recover  the  original  network.  Notice  that 
this  procedure  consists  of  adding  parasitic  capacitors  and 
inductors  at  appropriate  locations. 

In  order  to  give  a  transversalization  result  via  element 
perturbation,  let  us  first  define  a  (  ;  perturbation  of  A*. 
Let  M  be  a  C:  submanifold  of  R"  and  let  C:(  A/. 01”)  he 
the  set  of  all  C:  maps  from  M  into  R".  Let  A 6 (  :(  M . R ” ) 
and  consider 


-?l:(  /':«(■)) 


=  (,  :  Xf 


2  h  *  llt  <  My  . 

»  (41) 


(7  £  ('■(  M.U") 

if  A(  .r  )  -  (,{  x )  -  (JP)X  [d(i)x 
+  il (</->’),  (d:G  )r  <<(.r). 
for  a  I !  ,rt  M 


where  e(  x)  is  an  arbitrary  continuous  function  from  A/  into 
the  set  of  positive  numbers  and  d  p  and  d  (i  are  the 
second  derivatives.  These  sets  generate  the  strong  (  :  topol¬ 
ogy  for  (  ;(  M .  R " )  |6).  The  set  F.mb:(  M.U" )  of  all  (  ' 
embeddings  of  M  into  R”  is  open  with  respect  to  this 
topology  [6).  Let  Tl ■'(«„)  be  a  neighborhood  of  the  inclu¬ 
sion  map  such  that  all  elements  of  -?l '( i  „  )  are  embeddings. 
Then  a  C':  perturbation  A/  of  A/  is  defined  by  M  (Y(  A I ). 
where  "(j,, ).  The  following  is  our  first  transversah- 

z.ation  result  via  element  perturbation.  Although  the  proof 
is  similar  to  that  of  Theorem  3  of  (3|.  there  is  a  technical 
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difference  because  of  the  C2  perturbations  instead  of  (’* 
perturbations.  Proof  is  given  in  the  Appendix. 

In  the  sequel,  denotes  variables,  functions,  sets,  etc., 
associated  with  a  perturbed  network. 

Proposition  3:  Given  a  nonlinear  network  ‘.X  let  An  A' 
^  0  and  XSXK.  Suppose  that  A(  (resp.  A,)  is  locally 
voltage  (resp.  current)  controlled.  Then  there  is  a  perturba¬ 
tion  \R  of  A  H  arbitrarily  close  to  \R  in  the  strong  C2 
topology  such  that  A  fU'/0  and  A  ft  A,  where  A  is 
defined  by  (6)  in  which  A R  is  replaced  bv  Ay,. 

Remark:  Recall  Corollary  1  where  A,  (resp.  A ,  )  is 
locally  voltage  (resp.  current)  controlled  and  observe  that 
depends  only  on  ( vR,  ik ).  This  is  the  reason  why  one 
can  transversalize  A  and  K  by  perturbing  Ay,  only. 

The  next  result  gives  a  transversalization  procedure  via 
network  perturbation. 

Proposition  4:  Given  a  nonlinear  network  let  An  A 

0  and  AjJTA.  Let  '5  be  a  proper  tree  for  w  and  let  £.'  he 
its  associated  cotree.  Partition  V'  and  t  as  0  -  R  U  ('  and 
i  -  Rc  U  L.  respectively,  where  R.  C.  and  L  denote  resis¬ 
tors.  capacitors  and  inductors,  respectively.  Insert  an  arbi¬ 
trarily  small  linear  capacitor  in  parallel  with  each  branch 
of  R  and  insert  an  arbitrarily  small  linear  inductor  in 
series  with  each  branch  of  R  .  Then  the  perturbed  network 
\  satisfies  the  following  properties:  (i)  AnA'^fcP,  (ii) 
A  A . 

Proof:  Let  denote  the  branches  representing  the 
capacitors  added  in  parallel  with  R  and  let  A,  denote  the 
branches  representing  the  inductors  added  in  series  with 
R  .  Then  ~  ('U(',uR,;  is  a  proper  tree  for  'Tt  and 

t  /.  U  l. ,  u  R  is  its  associated  eotree.  Let 

*  (•'*„■  G  •  •'/,  ;  *Vr< ••'< ,)  | 

:  */?e"  *(  •  *V, 

4 


(43) 


we  have 

(«V  A  (47) 

where  C,  and  I. ,  are  the  capacitance  matrix  and  induc¬ 
tance  matrix,  respectively,  of  the  added  elements.  Since 
( rj(.  )  corresponds  to  open-circuiting  branches  of 

G,  and  short-circuiting  branches  of  /. ,  and  since  such  a 
situation  is  contained  in  A.  we  have  (47).  Next,  since  no 
resistors  are  added,  we  have 

A  "  {(  r.  i.q.4>)\(  v.  i.q.$)E  A  } .  (48) 

This  implies  that 

(  *'o  -  *()-<?«.-<»  n)e’\  (4t)> 

which  together  with  (47)  implies  (i). 

(ii)  In  order  to  prove  AffiA.  we  compute  *5)  ,f )  of  ( 17)  for 
Observe  that  fundamental  loop  matrix  B  for  associ¬ 
ated  with  the  tree  '7  assumes  the  following  form: 

*■>,  vt  r/ ,  ;  v/te  *< 

1 

1  «*, 

where  the  submatrices  in  (50)  correspond  to  those  of  H  for 
vt  (see  (23)).  The  sign  of  the  identitv  matrices  in  (50)  are 
chosen  merely  for  convenience  and  involves  no  loss  of 
generality.  Next,  notice  that 

/r~/k-  (*'«-*«)  (,V,k) 

».-,,/*•  °i  jx 

°lj*  D'Jh 

Substituting  these  and  (50)  into  ( 17)  we  have 


*< . 

1 

*/« 

Bkk 


(50) 


in 

K  /< 

*  C, 

D.ft 

1 

l\f, 

1 

be  the  variables  of  \  .  Let 


It  follows  from  ( 10)  and  ( 12)  that 


( *‘i,-  AHA  c 

(44) 

rank  |  OvJK 

;  D,kJr  ;  DikJx  ;  DiKJn  ln>..„ 

*'n  (*>••'/.  •  •'«  •  •'<  „ )  | 

(45) 

rank  ( D,  f,  : 

fi  l<  r,  .i/l  "l 

'll  (  ' K  ■  '  1  ,  ■  '#•'(,}  | 

We  first  claim  that  with 

rank  [»1(/;  . 

D<sfl  i(>,  <>)  "  1 

r„  (c„  r,  .0  :  r*  ■  f,  .  vR  ) 

for  all  (  r„.  iR  IfcA,.  (  r(  ,q)E.\,  and  (/, .  <f>)t  A 
lively.  It  follows  from  (51)  and  (48)  that 

'll  (  'it  •  '  1  '  R  'K 

(46) 

rank  ?(  x )  n  H  *•  n(  +  n,  +  n,  ♦  n, 

4n  (<)<,•£>>  )•*>  l4>o-L,iR  ) 


for  all  ie£.  where  //,  and  nt  are  the  number  of  capaci- 


II  l  I  IHANSAl  IIIINSOS  IlKUIIS  AMISISIIMS.  Mil  .  I  *S-’K.  Ml  MO  IWI 


■  W  '“4  • 
v  T  []- 

*  I _ 
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Fig,.  2.  A  nonlinear  circuit  with  Ai^A  (a)  The  circuit  diagram  (b)  Fig  F  A  nonlinear  circuit  whose  local  m>F  ability  depends  on  resisf«»i 


Resistor  constitutive  relations 


constitutive  relations  (a)  The  circuit  diagram  tb>  Resistor  constitutive 
relation  where  the  circuit  is  not  localls  solvable  tel  Resistor  constitu- 
live  relation  where  the  circuit  is  localls  solvable 


tors  and  inductors  added,  respectively.  It  follows  from 
Proposition  I  that  A  ft  K .  □ 

Example  2:  Consider  the  circuit  of  Fig.  2(a)  where  the 
resistor  constitutive  relations  are  given  in  Fig.  2(h)  with 
=  ).  k  =  1.2.  Then  5=  {/?,.  R:)  is  our  proper  tree 

and  BRR  =  BRI  -  Bt  /  -  0 Bt  r~[  1  1  ]. 


D,  fR  = 
»  «*  *  K 


n  fR  =  o 


matrix  for  X  is  given  by 


r«.  *r,  vi 


Therefore  (17)  is  given  by 
‘  -  />/*, 

tF(r.  /.<#>)- 


Now  for  the  value  i%  in  Fig.  2(h).  we  have  iR  -fR  {  vR[  ) 
-  Ip .( i>  v  )  and  (  DfRi  )t>  =  {D/r  ),k  =0.  Lei 

t  */»„  1  '  O  I’r  (r  )■  t/ n  -  "  lK  */.„  “  (R 

and  <{.„  =  Et  i  .  Then  {vR.vr.iR  .!,.<{>„)  ( r„.  <?>,,)£  2£ 

and  rank  ?(  r„.  /„.  <>„ ) :r  2  <  3  and  hence  \IBK.  Insert.  now. 
(  I,  and  C'p  as  in  Fig.  2(a).  Then  Proposition  4  tells  us  that 

A  KK. 

The  transversali/.ation  procedure  is  simplified  if  A  R  is 
locally  voltage  controlled  [3].  i.e..  (9)  holds  and 

rank  (  D1h  fR  )( (>  jj( ,  ~  n  R .  for  all  (  vR .  i H )  (£  A  R  n  i  R 

(52) 

or  locally  current  controlled,  i.e..  (9)  holds  and 

rank!  »,.„/*  )(1>  nR.  for  all  (  r„.  iR  )e  Ar  n  (  K. 

(53) 

Proposition  5:  Consider  the  situation  of  Proposition  4 
and  assume  that  Ar  is  locally  voltage  controlled.  Insert  a 
small  linear  capacitor  in  parallel  with  each  branch  of  R  . 
Then  the  perturbed  network  A  satisfies  the  follow  ing  prop¬ 
erties:  (i)  A  n  X" #  0  :  and  (ii)  A  iT  A. 

Proof  (i)  can  be  proved  jn  a  manner  similar  to  that  of 
Proposition  4.  It  is  clear  that  '  ■  CUC,  is  a  proper  tree  for 
A  and  i;  R  U  R  U  /.  is  its  associated  cotree,  w  here  C, 
represents  the  branches  of  the  capacitors  added,  lo  com¬ 
pute  ‘  i  (  v ).  observe  that  the  fundamental  loop 

' W ,■  Jcnolc  j  II  ■  0  matrix  hv 


B/o  Br  r 
B 1 1  Bj  r 


where  the  submatrices  are  those  of  .X .  Therefore. 
k'  •  A>f  fk  ®r,  /« 

A„,/r  ^  •  ^ik[  h  ;  fk  •  (r«.»R)  (*‘r-,r)- 

Substituting  these  into  ( 17)  we  have 


#(i) 


l>„f< 


By  using  (12).  (52)  and  elementary  operations,  one  can 
show  that  rank^(.r)  nR  •  n,  *  «(  •  n,  for  all  it  i'. 
where  is  the  number  of  capacitors  added.  Therefore. 
Proposition  I  implies  that  A  X K . 

A  dual  argument  shows  the  following: 

Proposition  6:  Under  the  same  setting  as  that  of  Pro¬ 
position  4.  assume  that  \R  is  locally  current  controlled. 
Insert  a  small  linear  inductor  in  series  with  each  branch  of 
R  .  Then  the  perturbed  network  X  has  the  following 
properties:  (i)  An K^=  o,  (n)  AiX K. 

IV.  Lot M  Sot  \. Min  in 

Recall  that  transversality  of  A  and  K  is  a  static  condition 
in  the  sense  that  it  has  nothing  to  do  with  the  dynamics  of 
X .  In  order  to  motivate  the  discussion  of  this  section  we 
first  consider  the  follow  ing  example. 

Example  3  Consider  the  circuit  of  Fig.  3(a)  where  Ar  is 
given  by  Fig.  3(b)  with  c#  aRt  iR  )•  If  we  choose  ■'  |('l  (o 

be  our  tree,  then  HK(  1.  ,  lK  *•  l>  ,  /r  IHr  ‘""-I 

rank ‘'4(  i  r.  Ir  )  rank)  I  '.  l)gR  |,  I  It  follows  from 
(’orollan  /  that  A  X  A  and  2:  is  a  perfectly  well-defined 
one-dimensional  submanifold.  Hie  dynamics,  however,  has 
points  where  it  is  not  well  defined,  lo  show  this  observe 
that  /r  serves  as  a  global  coordinate  for  1.  i  c . 

(  I  r.  I',  -  I  R.  I,  .  (/  )  (  Sr(  1  R  )■  V’r(  I  R  ).  1  R.  1  R.  3  Cr<  I  R  I). 

where  (  is  the  capacitance.  (Nonce  that  r,  cannot  serve  as 
a  coordinate.)  In  terms  of  this  coordinate,  the  dummies  is 
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given  by 

C(Ds*]"lf  =  "'«-  (54) 

Since  (Dg„), s  -(£>£*),,,  =  0,  differential  equation  (54) 
is  undefined  at  iR  =  /R  and  iR  =  /„  .  Then  one  might  like 
to  choose  another  coordinate  and  check  if  the  dynamics  is 
well  defined  everywhere  in  terms  of  it.  If  it  fails,  then  one 
may  try  to  choose  another  coordinate  and  repeat  the  same 
procedure.  The  problem,  however,  is  that  there  are.  in 
general,  infinitely  many  coordinates.  Therefore,  well- 
definedness  of  network  dynamics  should  be  defined  in  a 
coordinate-free  manner  and  there  should  be  a  method  of 
checking  that  porpertv  in  a  coordinate  independent  manner. 
That  is  exactly  the  problem  of  local  solvability  discussed  in 
this  section. 

We  will  first  show  how  the  dynamics  of  a  network  is 
described  in  a  coordinate-free  manner.  Let  it':  R:/’'"'  *"< 
— »R''>  *">  be  the  projection  map  defined  by 

ir '(t\i.?. £)  =  (</,<>)  (55) 

and  let 

w=  ir'oi  (56) 

where  i  is  defined  by  (21).  This  map  is  the  same  as  the 
restriction  ir'|22  of  it'  to  2.  Consider  the  following  symmet¬ 
ric  two-tensor  G  on  R ”<  *”>: 


two-tensors  on  a  general  manifold  since  manifolds  gener¬ 
ally  are  nonlinear. 

Let  ir*  be  the  induced  map  |6]  of  it.  Then  w*  pulls  G 
back  to  2  by  the  following  formula: 

(w*G)x(£,.£:)=  Gv(J{di,)J,.{dv)xi:)  m 

where  &  TX^,.  Similarly,  let  i*  be  the  induced  map  of 
i.  Then 

to-  i*rj  (60) 

is  a  one-form  on  2:  defined  by 

«,(£)-  n.,x ),!;)■  (61) 

A  coordinate-free  description  of  the  dynamics  is  given  by 
the  following: 

Proposition  7:  Let  —  be  an  (  n(  +nt  (-dimensional  C' 
submanifold.  Then  the  vector  field  A0  describing  the  dy¬ 
namics  of  network  satisfies  the  following: 

(ir*G  ),(*,.  4)  =  to,U).  forall4ertv.  (62) 

Remark:  An  easy  wav  of  understanding  (62)  is  the  fol¬ 
lowing:  Let  A  be  a  symmetric  n v  n  nonsingular  matrix. 
Then  (Ax.  y)  is  a  symmetric  bilinear  function  on  R”  -  R”. 
where  (•■•)  is  the  inner  product,  i.e..  (Ax.  y)  ( x.  A  'y ). 
and  (/4.r.  •)  and  (/j-.y)  are  linear.  For  a  vector  to£R". 
the  formula 


nl  « / 

G-  2  t/t/*®^  _  2  d<t>k®d< Pk  (57) 

k  \  k  - 1 

and  the  following  one-form  onR:'”"1  * : 

M,  H  i 

»i~  2  'i\d<ik  -  2  (58> 

it  i  k  i 


Remark:  A  simple  explanation  of  one-forms  is  given  in 
1 3 J.  A  symmetric  two-tensor  G  on  R:  is  a  collection  of 
functions:  R:x  R:-*R  given  at  each  point  (v,..v:)€R: 
by 

Gi  *,.«,)=  2  /m„(v,..v;)d.v„,(S'd.v„ 

m.n  \ 

where  are  real-valued  functions.  and 


dx ,  ®  dx ,  = 

1 

.  t/.Y,®t/v,= 

•  1 

dx2®dx,  = 

1 

.  dx2®d.x  2  = 

Consequently. 

1 

G|  > , . , ,  i(l  I  0lr. 

(1  0J7)  =/,,(., 

1  •  ) 

where  we  look  at  a  matrix 
satisfying 


b 

d 


as  a  map:  R:  xR: 


—  R 


a 

( 


a 

h 

1  v:  l 

l) 

c 

d 

1  I 

Therefore.  G  can  be  thought  of  as  the  matrix-valued  func¬ 
tion  |  fma].  One  needs  to  be  careful,  however,  in  defining 


(  Ax.  y)  =  (<o.  y).  for  ally  PR” 

uniquely  defines  the  vector  x  A  'to.  If  the  network  is 
reciprocal  and  if  P  denotes  the  mixed  potential,  then  to  of 
(62)  is  given  by  t c  =  dP.  differential  of  P.  If  tq.  <>)  serves  as 
a  global  coordinate  system  for  12.  then 


-I,  n, 

«,  =  2  /V.(<?-<M‘7‘/t  2  /  /.(</•<>) 

A  1  A  I 

where  F(  and  F,  are  delermined  by 

(»<  •  -«•'/ ) 

-  ( ^;, (?•£)•■  .<?•<!>)•  Fi  ,<  ■<!•)■■■■ -i) ,  (q$))- 

Therefore.  (62)  is  reduced  to 

^ (?•£).  k  >  •■••«< 

~j*-  ~ ■/■).(*■£)•  k  )•  ••«/• 

Proof  of  Proposition  Let 

r  (63) 


A  vector  fickl  V  on  manifold  2  is  ,i  function  such  ih.it  the  value  V,  .11 
.rf  2  belong  to  l\  2 .  the  tangent  space  of  2  .it  .1  T  he  set.  lot  fiekl  \ 
n.ituralh  generates  a  floss  .v| f )  such  that  </.vtn  Ji  \xl., 


L, 
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Then,  by  definition,  the  left-hand  side  of  (62)  reads 
{ m*G  )  x(  Xx .  i)  -  G„,  x\(dm)Xx .< dm )  xi ) 


-  2  dqk^dqk 


2  d<j>k®d$k 

k  1 


j  dq  <Af>  \  \ 

J  <"  '  d>  I 

y  <j<h,  _  y  «/$*, 
^  dt  ^  L  '■  ^ 


dt 


(64) 


The  right-hand  side  of  (62)  is  given  by 

n. 

2  '< .  %  - 

A  l  A 

W(  «/  \ 

2  '(  .«'</*  -  2  ‘y/Aa  «</*),*) 

A  1  At  / 

n(  n, 

=  2  A. f,(-  2  IT AV  (65) 

*  i  it  i 

Equations  (64)  and  (65)  together  with  (3)  and  (5)  implv 
(62).  □ 

Proposition  7  and  Example  3  naturally  lead  to  the  follow¬ 
ing  definition: 

Definition  J:  Given  a  nonlinear  network  /U  assume  that 
2!  is  an  (n<  +  nt  (-dimensional  C2  submanifold.  Then  ,Yt  is 
said  to  be  locally  solvable  if  (62)  uniquely  defines  a  C1 
vector  Xx  S  7X2  at  each  point  jc€2. 

Remark:  Local  solvability  defined  above  is  a  coordinate- 
free  version  of  the  one  in  Chua  and  Wang  (7|.  If  A(  (resp. 
A,  )  is  globally  voltage  (resp.  current)  controlled,  this  defi¬ 
nition  coincides  with  regularity  of  Smale  (8). 

We  are  now  ready  to  stale  the  main  result  of  this  paper. 

Theorem  I:  Given  a  nonlinear  network  vU  suppose  that 
X  is  an  («<  +«(  (-dimensional  C2  submanifold.  Then  \Yt  is 
locally  solvable  if  and  only  if.  at  each  point  xEZ,  (q.<j>) 
serves  as  a  local  coordinate  system  for  22. 

For  proof  we  need  a  lemma.  Recall  ( it* O’  )*(•••)  defined 
by  (59)  is  a  bilinear  function  on  7,22  X  7,22. 

Lemma  I:  Suppose  that  2  is  an  (n(  +n,  (-dimensional 
C;  submanifold.  Then  is  locally  solvable  if  and  only  if 
at  each  point  jre2.(ir*G'  )x  is  nonsingular,  i.e.. 

(w*G’)j(4,.{:)  =  0,  for  all  €7^2  implies  £,  =0. 

(66) 


<t*n)jr(€)=I 


((*),{) 


Proof:  We  look  at  (ir'G)^-.  • )  in  a  slightly  different 
manner.  Consider  the  map  Jx  defined  by 

G)x({,.-).  (67) 

To  each  £,.  the  map  Jx  assigns  the  linear  functional' 
(w*G)x(£,.  •)  on  7,2,  A  linear  functional  on  Tx 2  belongs 
to  its  dual  7/2.  This  means  that  Jx  ( ir*G ),( • ,  )  maps 
Tx 2  into  7/2.  It  is  clear  that  (66)  implies  that  Jx  is  an 
isomorphism  and.  therefore,  it  is  invertible.  It  follows  from 


'A  linear  functional  is  a  real  valued  linear  function 


(62)  that  a  vector  field  Xx  is  uniquely  determined  by 

Xx^(Jx)  (68) 


In  order  to  show  that  X  is  C'1.  recall  definition  (59)  of  m*(i. 
Since  2  is  C2.  ihe  map  dit  is  C1.  Therefore.  (  Jx )  1  is  C1 
Similarly  to  is  C1.  This  implies  that  X  determined  by  (68)  is 
C1.  Conversely,  if  J,  is  not  an  isomorphism  (62)  cannot 
determine  a  unique  vector  field. 

Remark:  In  Example  3.  in  terms  of  the  coordinate  iH.  we 
have  (ir*G  ),  =  C(  Dv,R ),  diR<S>diR  w  hich  becomes  singular 
when  ( l)yH  ),a  =  (  Dyk  )*  0. 

Proof  of  Theorem  I:  Recall  definition  (59)  of  ir*G. 
Since  G'„(t)  defined  by  (57)  is  always  nonsingular  in  the 
sense  of  Lemma  /.  we  see  that  (  m*G  ),  is  nonsingular  if  and 
only  if  the  following  map  is  an  isomorphism: 

(dm)x:  7,2  -  7„(„R"'  -<  (69) 


i.e.,  it  is  a  local  diffeomorphism  al  ar.  But  this  precisely 
means  that  serves  as  a  local  coordinate  system  for  2 

at  x.  Zi 

Remarks:  l )  Because  of  its  coordinate-free  property  . 
Theorem  I  is  of  fundamental  importance.  It  says  that  if  all 
the  variables  of  a  network  are  expressible  in  terms  of 
capacitor  charges  and  inductor  fluxes,  at  least  locally,  then 
the  network  is  locally  solvable,  i.e..  the  dynamics  is  well 
defined.  Conversely,  if  the  network  is  locally  solvable,  then 
necessarily  determines  all  the  variables  in  the  net¬ 
work.  Another  important  interpretation  of  Theorem  /  is 
that  if  (62)  uniquely  defines  a  C1  vector  field  with  respect 
to  one  coordinate  system,  then  it  defines  a  unique  C' 
vector  field  with  respect  to  capacitor  charges  and  inductor 
fluxes  also.  Conversely,  if  (62)  fails  to  specify  a  unique  C1 
vector  field  with  respect  to  j>).  then  there  is  no  choice  of 
variables  in  the  network  in  terms  of  which  (62)  specifies  a 
unique  C’  vector  field.  These  observations  show  that 
capacitor  charges  and  inductor  fluxes  are  the  fundamental 
quantities  in  describing  the  dynamics  of  a  network. 

2)  Let  us  explain  why  A  must  be  C2  in  order  to  define  a 
C1  vector  field  by  using  a  simple  example.  Consider  the 
circuit  of  Fig.  3(a)  where  A „  is  given  bv  Fig.  3(c).  Assume 
that  g„  is  a  global  C1  diffeomorphism.  Therefore  iR 
hR(vR).  where  hH  =  jt,R  1  and  hR  is  also  a  global  (  1  diffeo¬ 
morphism.  The  sets  A  and  2  are  C’  submanifolds.  Capaci¬ 
tor  voltage  r,  serves  as  a  global  coordinate  for  2  and  the 
dynamics  is  given  by 


di\  ^  hR(  v,  ) 
dt  C 


(70) 


The  right-hand  side  is  C1.  Now  it  is  clear  that  iR  is  another 
global  coordinate  for  2  and  the  dynamics  is  given  by 


djjt  _  *  w 

dt  C(T>g*),/ 


(71) 


Since  gR  is  C1.  the  right-hand  side  is  only  C".  This  gives 
rise  to  a  problem  because  a  C"  vector  field  cannot  guaran¬ 
tee  uniqueness  of  solutions.  If  we  assume,  however,  that  gw 
is  a  C2  global  diffeomorphism.  then  the  right-hand  side  of 
(70)  and  (71)  is  at  least  C*.  Therefore.  C'-ness  does  not 
depend  on  the  choices  of  coordinates.  More  generally,  let  X 
be  a  vector  field  on  2  and  let  (L/n2,^)  be  a  local  chart  at 


je. 


-  ■  ^  — 
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x6l  Then  a  natural  coordinate  representation  is 

(*<*»*„  =  <**>,*.■ 

If  is  another  chart,  then  for  xei'fU'ni’  one 

has 

Therefore,  if  we  want  Af(i£)  to  he  C'  independent  of  the 
choice  of  coordinates,  we  must  require  the  change  of 
coordinates  ')  to  be  Cr.  This  requires  tj>°  1  to 

he  at  least  Cr‘ '.  But  this  is  exactly  the  condition  required 
for  2  to  he  at  least  C * Therefore,  a  C  vector  field  can 
possibly  he  well-defined  only  on  Cs  manifolds  with  .v  >  r . 

3)  Observe  that  C'-ness  of  vector  field  is  required  in 
order  to  guarantee  uniqueness  of  flows  because  a  C°  vector 
field,  in  general,  does  not  suffice  for  generating  a  unique 
flow.  Roska  (9]  obtained  several  uniqueness  results  in  terms 
of  the  network  topology  and  resistor  constitutive  relations. 

4)  If  dim^!^t«( ■  +  n, .  then  (dir)x  of  (69)  is  always 
singular  and  (62)  cannot  determine  a  unique  vector  field. 
Therefore,  transversality  of  A  and  K  is  one  of  the  im¬ 
portant  conditions  for  local  solvability. 

The  following  shows  that  the  results  of  (5|.  [8|,  and  ( 10] 
are  a  special  case  of  Theorem  I . 

Corollary  4:  Suppose  that  At  (resp.  A,  )  is  represented 
by  q- R<(t-(  )  (resp.  <i>~g,ii, ))  and  suppose  that  (Dg,  ),. 
(resp.  (Dg,  )l( )  is  nonsingular.  Let  2  be  an  (n(  +nr)- 
dimensional  C2  submanifold.  Then  y't  is  locally  solvable  if 
and  only  if  at  each  point  jr  G  22  (e,  ./,  )  serves  as  a  local 
coordinate  system  for  2. 

Proof:  If  (Dgc)Vi  (resp.  (Or,),,)  is  nonsingular.  g, 
(resp.  g, )  is  a  Ureal  difleomorphism.  Therefore,  (£<•.«',) 
serves  as  a  Ureal  coordinate  system  for  2]  if  and  only  if 
(q.  <[> )  serves  as  a  local  coordinate  system  for  2.  □ 

Remark:  In  [5],  ]S).  and  [10],  ( Dgt  )t.  and  ( Dg, )(/  are 
symmetric  and  positive  definite.  Therefore,  they  are  non¬ 
singular. 

The  following  is  an  example  of  a  locally  solvable  circuit 
whose  capacitor  is  not  voltage  controlled. 

Example  4:  Consider  the  circuit  of  Fig.  1(a)  where  A{  is 
given  by  Fig.  1(c)  and  \R  is  given  by  tR  -  gR(  vR).  Capaci¬ 
tor  charge  q  is  a  global  coordinate  for  2  and  the  dynamics 
is  described  by  q-gR(~gc{q)).  Clearly,  this  circuit  is 
locally  solvable  but  the  dynamics  cannot  be  described  in 
terms  of  v( . 

We  will  next  show  that  if  the  assumptions  of  Corollary  4 
are  satisfied,  then  (62)  is  reduced  to  a  formula  in  [10].  To 
this  end  let 

2f  =  A tnKt  (72) 


where 

A*  {(f.i)eR^|(tvt,)€A#j  (73) 


K*  =  ((r.i')€R:*|(r.i)  satisfies  KVL,  KCL}.  (74) 

Let  A,  and  A,  be  characterized  as  in  Corollary  4  and 
define 

n(  n  i 

2  Cmn(  t>,  )  dv(  _  ®  dx\  _  -  2  Lmn(i,)dt,®di^ 

ni.ti  I  m .  n  I 

where 

|UOl  [/-m„(c  )]  -  (OR/  )„.  (75) 

Let  — *R'''  '  be  the  projection  map  defined  by 


wf(  r. « )  -  ( rt  ) 


(76) 


and  let  t*:  2!1  — •R"''  be  the  inclusion  map.  Finally,  let 

V  r  2  +  </  2  '(  /(. 

A  I  \  k  I 

Then  the  vector  field  .V*  which  describes  the  dvnamics  is 
given  by  [  10] 


(wt*Ct),r.o(Ao..v^)  =  <,,(r).foralUte7;i,112:1 


(77) 

where  uf  =  t1*!)1.  Let  F.  —>2;  be  the  (global)  dilfeomor- 
phism  defined  by 

F(r.i)  =  («•.«, *<(»•<  ).gi(i,  )).  (78) 

Proposition  8:  Suppose  that  the  assumptions  of  Corollary 
4  are  satisfied.  If  '.A  is  locally  solvable,  then  (62)  is  reduced 
to  (77).  i.e., 

K  (79) 

Proof:  Let  Xf{r  be  the  vector  field  determined  by 

(62).  It  follows  from  (78)  that 

Xnr  .o=(dF\rJ)Xt\-.o  m 


for  some  AT,*, = 

x:={x:„ 

Then 


( X f,  X*)eT{r Let  us  write 

-  •  Xv, )  ■  X*  r~  ( A',*  •  Af,* .  A  7  ) . 


( dF  )(r,,V,,o  =  {K'  X* .( %<),.,  Xl  .( Dg,  ),.  AT,*  ) . 

We  will  show  that  X,ri)  is  the  same  as  the  one  determined 
from  (77).  To  this  end  note  that  for  £G  Tx1,  there  is  a 
G  T{r  ,-,2*  such  that 


t=l4F\.'„V  =[&tUDgc)r,£,  -{Dg,  ).(l  ).  (81) 


We  substitute  (80)  and  (81)  into  (62).  Then  the  left-hand 
side  reads 


k  i 


K  I 


2  dqk®dqk~  2 

2  (  2  i(  2  2  (  2  2 

m  "  I  \  n  --  I 


m  I  *  n 


n  I 


nit  }  K  \  NS  U  J  U  »Ss  UN  t  1K‘  l  M  s  \M  >  M  Si  I  MV  ^  ‘  »l  <  -wZN. 


M\i  ;  l'.s  1 


where  we  used  (75).  The  right-hand  side  Kirns  out  to  he 


m  1  '  n  I 


£  v.  £  £  /• . (<i  h!,  ■ 


ni  I  '  u  \ 


Si  nee  f  ( d(  rr  since  |i//')(t  rl  is  nonsingular  and 

since  is  locally  solvable,  i.e..  {dtt)tu.  h  is  nonsingular, 
there  is  a  | 1  such  that 

«/ 

£  <;,„(«•«  )*,*.  A .  £  )il  0  (84) 

n  I  n  1 

w  here  Snw  I  if  m  n  and  Snw  0  if  m  a  n 
Substitution  of  (84)  into  (82)  and  (83)  gives 

>/, 

£  )-v,.  i»5) 

n  1 

Similarly,  there  is  another  vector  >ueh  that 

n  i 

£<;,„<«<  . (h)'K,  s”--  <86> 

n  i  //  i 

Substituting  (86)  into  (82)  and  (83)  we  have 

£  )  V„  (87) 

*i  I 

The  vector  field  ,Vtl. ,,  determined  by  (85)  and  (87)  is 
exactly  the  same  as  the  one  obtained  by  (77).  2 

We  will  next  give  a  simple  formula  for  checking  local 
solvability  of  \ .  Let  /„.  /  .  and  /,  be  as  in  (9)  ( 12). 

Proposition  v.  Let  -  be  an  <«,  *  «,  )-dime»siona!  C  : 
submanifold.  Pick  a  proper  tree  ■'»  and  let  c  be  its  associ¬ 
ated  eotree.  Then  is  locally  solvable  if  and  only  if 

det9C(.v )=*=(.).  for  all  .t£22  (88) 

where 

[»,,/«  :  Di(fR  4  1 


where  B.  Q.  etc.,  are  as  in  1 18).  Since  22i  (  (i  '(0).  we 
have  (31 

7,1'  lm{P4  ’I*,,  Ker(/W),  (93) 

It  follows  from  (91)  and  (93)  that  the  matrix  of  (90)  i-. 
nonsingular  if  and  only  if 

Ker(  />«■’),(''  Ker(  D(i  ) ,  0  (44) 

which  is  equivalent  to 


2 />  +  u,  •  n /  . 


Computing  the  matrix  of  (95).  one  can  show  that  it  has 
rank  Zb  *  n(  -*  n,  if  and  only  if  the  following  matrix  has 
rank  Zb: 


t'r 

tv 

#L  / 

1 

B 

(86) 

B'  1 

l\  ,  fx 

/>,./*  />, .  (h 

(87) 

:  f,c  h 

By  elementary  operations,  one  can  show  that  this  matrix 
has  rank  2 b  if  and  only  if  (88)  holds. 

fins  result  has  an  interesting  consequence.  Let  n,  2.'  - 
IR;"‘  (resp.  it,:  22  — R:'’'  )  be  the  projection  map  defined 
b\  it,  U)  (c,  .q).  (resp.  ir,t.r)‘  («',  .$)) 

Corollary  5.  If  is  locallx  solv  able,  then 


J*  (!>,,/<  )„.  v, 


Jet  [OJ, 


for  all  (  r,  .*/)•.  it,  ( 22 )  (9') 
for  all  ((,.<>  |  c:ir(  ( 1' ).  ( 98 1 


».,/«  (  '>,,/*  >« 
K  /< 


Proof :  Let  (^,220  L' )  be  a  local  chart  for  22  at  .v  Then 
( dv ),  is  an  isomorphism  if  and  only  if  (  0(  tr<’  1 1)^, ,,  is  a 
nonsingular  matrix.  Since  ^  1  we  have 

(»(* '))*«.  h !>*'),(*),( ©+  m 

Since  (d t),  is  a  linear  inclusion  map.  the  matrix  of  (90)  is 
nonsingular  if  and  only  if 

Ker(0i»')tnim|/>*  |0).  (*M ) 

Let  (i:  (.  -R;"  be  defined  bv 


Qi 

f»  (  AT  )  /*(  •>  •  *R  ) 

/<  ( *'(  R  ) 

/,('/•£) 


Proof:  Since  (7?,.  /  ),,.  and  ( />,  /,  ),,  (  are  square 
matrices.  (88)  forces  (97)  and  (98)  to  hold. 

Remark:  The  above  result  savs  that  if  \  i s  localh 
solvable,  then  capacitors  must  be  localh  choree  controlled  on 
ir,(22)  anil  inductors  must  be  localh  flu  \  tonrrol/cd  on 
11,(22).  In  other  words,  if  capacitors  are  not  locally  charge 
controlled  at  some  point  ( r,  .</,,)Lir,  (22)  or  inductors  are 
not  locally  flux  controlled  at  some  point  (i,  .  <>„  K:  it,  i  22 ). 
then  there  is  no  choice  of  local  coordinate  s' stem  m  terms 
of  which  the  network  is  locally  solvable.  Ihts  seems  to 
explain  whv  we  do  not  find  capacitors  (resp.  inductors i  of 
practical  importance  which  are  not  locally  charge  (resp. 
flux)  controlled. 

Importance  of  local  charge  (resp.  flux)  eontrolledness  is 
further  emphasized  by  the  following: 

Corollary  b:  Suppose  that  A,  (resp.  A,  |  is  described  bv 
q  ~  g(  (  r,  )  ( resp.  <>  r,(  i,  »  and  suppose  that  t  r,  .  i,  I  serv es 
as  a  global  coordinate  system  for  22  If  det(  1)%,  )r,  0  for 

some  r,  tR"'  or  del  ( /)j»,  )(  0  for  some  i,  i  R”  .  then 

A  is  not  IvK'allv  solvable. 

Proof:  If  (r(  .  i,  )  serves  as  a  global  coordinate  sxstent 
for  22.  then  ( v(  .  g,  ( «;  ( )  t.  -jt,  <  22  >  for  all  r,  t  R"  aiivl 
(/,.£,(/, ))  £11,(22)  for  all  »,  6R'1'.  Therefore.  ( 'ort diary  5 
implies  the  result. 
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Observe  that  the  above  result  says  that  even  it  ( r(  . ) 
serves  as  a  global  coordinate  system  tor  2A  a  network  mat 
not  be  locally  solvable.  The  following  example  shows  the 
ease  in  point 

Example  5.  Consider  the  circuit  of  f  ig.  1(a)  where  A,  is 
given  bv  Tig.  1(c).  As  was  shown  m  Txample  I.  A  T  A 
always  holds  for  am  \A,  as  long  as  it  is  a  C*  submanifold 
Suppose  that  \K  is  described  bs  /A.  gA(rA. ).  I  hen  i, 
serves  as  a  global  coordinate  for  2A  Since  there  are  points 
where  t  /)tr,  ),  t>.  this  circuit  is  not  locallv  solvable. 

If  we  know  that  A,  tresp.  A ,  )  is  locallv  charge  (resp. 
flux)  contri'lled.  then  the  formula  in  Proposition  '>  is  sim¬ 
plified  as  follows 

Corollan  ~ .  In  the  same  setting  as  m  Proposition  v 
assume  that  \,  (resp  \,  i  is  locallv  charge  tresp.  flux) 
controlled.  Then  A  is  locallv  solvable  if  and  onlx  if 

del  0CA.|  rA .  i"  )  *  0.  for  all  (  rA. .  /  A. )  •  irA.(  —  I  I1)4) 

where  irA.  is  defined  In  (22)  and 

rn  •  'k  ) 

[ /h  (  >KJh  )bkh  :  />,.  /,  •  ( i\JR  )^4,,  , 

(  100) 


Proof:  If  the  above  hypothesis  is  satisfied.  (  /),  /  ),, 
and  (/>,/,  tire  nonsingular.  Then  one  can  show  that 

(00)  is  equivalent  to  (88). 

Remark  In  [5|.  [8].  and  |10).  A,  (resp.  A,  )  is  repre¬ 
sented  bv  q  g(  (t-  )  (resp.  <5  gpi,  ))  and  ( /)g(  ),  (resp. 
I  D g,  )„  )  is  positive  definite,  Therefore.  A,  (resp.  A,  )  is 
locally  charge  (resp.  flux)  controlled. 

Example  b:  Consider  the  circuit  of  Txample  V  Since 
5A(r R.iR)  (A)##),  .  it  fails  to  have  rank  1  at  / A,  t A. 

and  iK  i K  and.  therefore,  this  circuit  is  not  locallv  solva¬ 
ble. 

Example  Consider  the  circuit  of  Txample  2.  where 
\R  is  given  in  Tig.  4(a)  with  iK  XrU'r  ^  1.2.  Since 


vr ■ 'h  ) 


1 -ji  4  A  iionlmc.u  circuii  which  in  noi  1«h.jI!\  M>l\»thlc  uti  RcvMi 
ki'iiNtitutiNC  relations  i  b)  Protection  of  1  onto  the  t  r  ,  . .  ,  i  n|\uc 


where  o.  />.  y.  and  6  are  as  in  (24).  In  particular,  if  \A.  is 
sloballv  volume  controlled.  then 


0CA.(  rA.  i  A.  ) 


I  DE 


B, 


1 

Bit, 


(102) 


and  if  \A.  is  ghiballx  current  controlled,  then 


lR  •  ‘r  I 


(  DC  I 


Bin 


(  10  .' I 


Recall  (24).  (35)  and  the  notation  used  m  Proposition  3. 
Proposition  It i.  l  et  Aa.  A,  .  and  A ;  be  globally  param¬ 
etrized  and  let  2!  be  an  ( n,  •  n ,  (-dimensional  (  :  submani¬ 
fold.  Pick  a  proper  tree  and  let  >.  be  its  associated  cotree. 
Then  A  is  locally  solvable  if  and  only  if  for  each  pCR' 
w  ith  .r(  p.  /,  .  v,  )-l  21. 

detDC*(p|x(l  (1041 

w  here 
Ot*(p) 


Bvh 

*  Bhh(  ) 

Br  ,  (  «.<  ) 

Bin, 

BiAOiRj 

Dt, 

«/A.(  Di,  I 

/><> 

(1051 

Proof:  Substitute  (34)  and  (35)  into  K.VL  and  kCT: 


and  since  PgK  and  OgH  never  vanish  simultaneously,  rank 
&i(  t'R.  iR )  ~  2.  Consequently  AT  A  and  2]  is  a  one- 
dimensional  submanifold.  Since 

1k-  ) 

there  are  points  w  here  del  9C^(  vx.  iK )  -  0.  Therefore,  the 
circuit  is  not  locally  solvable.  If  we  use  Corollary  4.  we  can 
see  this  more  clearly.  Consider  the  projection  ‘A  of  21  onto 
the  (  r, .  Tspace  given  in  Tig.  4(b).  If  we  further  project  'A 
onto  the  r,  -axis,  we  see  that  i,  cannot  be  a  local  coordinate 
where  the  curve  intersects  itself.  Therefore.  .A  is  not  locally 
solvable. 

Corollary  H:  Let  A  R  admit  a  generalized  port  coordinate 
representation  and  let  A,  (resp.  A,)  be  locally  charge 
(resp.  flux)  controlled.  Then 


1  ;  BhsHh^ 

1  B,  hBu 


<■«,<  Pr  I 


•'rJPr  ) 
•  '<  <  P<  ) 


(10b) 


Bit,  Bin  i  I 

BL  «/, ;  ■ 


l  et  us  write  (l()b)  and  (107)  as 


IrSPr  * 
'/(Pi  I 
I  IrAPr) 

>i 


//(  p.  i,  .  r,  )  0. 


0 


(107) 


(108) 


I  hen  2  is  diffeomorphic  to  //  '(0)  Bv  a  similar  argument 

to  that  of  the  proof  of  Proposition  *A  one  sees  that  (4(1)  is 


%„(«>.••*)  [(«,  (PE)y2) -(«,  (DE)yx)BHH  \  (ft  <  W  I*,)  M  ft  (  OF  )S:  )/»«'*]„  ^  , 


(1011 
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♦  *f-i  *  -  *  *<-  1  - 


r- — — 

r — (y 

f1' 

•q 

■c. 

l-f-h 

£~ > 

-  *R,  *  * 

Fig  5.  A  non/incar  circuit  which  is  iocailv  soluble,  vet  -  is  not 
diffcomorphic  to  Rn 


nonsingular  if  and  only  if 


rank 


DH 

Dv' 


=  b  +  n(  +  n  r . 


(109) 


I!  /(|  v).  v.  r  )1i  — >  x  as  ii(  y.  j  )-  — *  x.  Hence  F  is  proper  [6], 
Consequently  it  is  an  embedding  |6).  Define 

A*  ■  F(  S'  ,<R;) 


Computing  the  matrix  of  (109)  and  using  elementary  oper¬ 
ations,  one  can  show  that  ( 109)  is  equivalent  to  ( 104).  □ 

Remark:  The  above  result  holds  even  when  \R.  A(  and 
A,  are  locally  parametrized. 

Now,  observe  that  OC(jt)  of  (89)  is  a  submatrix  of  €(.  ' 
defined  by  (17).  This  implies  the  following: 

Proposition  11:  If  (88)  holds,  then  AffiA  and  X  is  lo¬ 
cally  solvable. 

Remark:  A  similar  result  holds  for  Proposition  10  and 
Proposition  2.  Observe  that  while  Proposition  9  assumes 
that  ^  is  an  (  nc+n,  (-dimensional  C:  submanifold.  Pro¬ 
position  11  does  not. 

In  many  practical  networks,  it  is  a  global  diffeomor- 
phism.  i.e..  all  variables  in  the  network  can  be  globally 
expressed  as  a  function  of  (q.  tj>)  and  hence  £  is  globally 
diffcomorphic  to  R"1  * .  Of  course  '.X  is  locally  solvable. 
In  Example  7.  "X  is  not  locally  solvable  and  £  is  not 
diffcomorphic  to  R.  A  question  arises;  Are  there  networks 
such  that  2  is  a  submanifold  not  diffeomorphic  to  R'1'  ’ 
yet  they  are  locally  solvable?  The  answer  is  affirmative  as 
the  following  example  shows. 

Example  8:  Consider  the  map  F:  RJ  —  R*  defined  by 
F(x.  y,  r )  ~  (e'eos.v.  e'  sin  v.  y. cos  v.  sin  v).  For  v.  ,v*€ 
R.  define  the  equivalence  relation  x~~x'  bv  ,v-.v'-2 kit 
where  k  is  an  integer.  Clearly,  then,  the  quotient  space  of  R 
with  respect  to  this  equivalence  relation  can  be  regarded  as 
the  unit  circle  S1  in  R2;  R/~=S'.  Let  (v)  denote  the 
equivalence  class.  Then  F  naturally  induces  the  map  F: 
S'  XR:-,R‘  by 

F([.x],.v,r)  -  F(x.y.z).  (110) 

Since 


e'(  — sin.v)  e'eos.v 
e'(cosjv)  e'sinx 

1 

-sin  x  cos.v 


-111  <|.  v.i) 


and  since 


det 


ev(— sin.v)  e'eos.v 
ev(cos.v)  e'sinar 


eiv^0 


iR  e'eos.v.  iR  e '  sin  v,  iR  : 
i  R  i.  tH  com.  ;a  sin.v.  (Ill) 

This  is  a  parametric  representation  of  \  R.  Consider  the 
circuit  of  Fig.  5  where  \  R  is  described  by  (111).  It  follows 
from  the  above  argument  that  \  K  is  a  three-dimensional 
submanifold  diffeomorphic  to  A'1  >  R  :.  It  is  clear  that  )£  is 
diffeomorphic  to  \R  and  therefore  diffeomorphic  to  .V1  ■ 
R;.  Notice  that  pR  —  (  v.  y  . : )  always  serves  as  a  local 
coordinate  system  for  A R  (not  a  global  coordinate  system, 
however).  As  we  remarked  earlier.  Proposition  10  holds 
even  when  \R  is  locally  parametrized.  The  matrix  of  (105) 
is  given  by 

9C»(p> 


—  e '  sin  v  e '  cos  v 
e'eos.v  e' sin.v 


Ct 


1 


p 


Since  detOC*(p)  «**’  =X=0.  X  is  locally  solvable.  Conse¬ 
quently.  the  dynamics  of  X  is  perfectly  well  defined  on  2. 
yet  there  is  no  global  coordinate  system  in  terms  of  which 
the  dynamics  admits  a  global  state  equation  because  2! 

S'  xR:  ^R'. 

Next,  we  will  give  two  more  examples  that  are  of  inter¬ 
est. 

Example  V.  This  example  shows  that  there  is  a  nontriv¬ 
ial  locally  solvable  circuit  whose  inductor  is  locally  flux 
controlled,  but  not  locally  current  controlled.  Consider  the 
circuit  of  Fig.  6(a)  which  consists  of  a  I -ft  linear  resistor 
and  a  Josephson  Junction  device  characterized  by  /,  - 
A  ,  sin  k:<t>.  where  A,  and  A ,  are  constants  (Fig.  6(b)).  This  is 
a  flux  controlled  inductor  which  is  not  locally  current 
controlled.  One  can  easily  show  that  transversalitv  and 
local  solvability  are  satisfied. 

Example  10:  Here  we  will  illustrate  a  power  of  geomet¬ 
ric  approach  using  an  interesting  example  of  Gbcknar  [)  1). 
Consider  the  circuit  of  Fig.  7(a)  where  the  resistor  is  linear 
1-0  and  A,  is  characterized  by 


we  have  rank  (dF)i |t), =  3  for  all  ((.v],>,;)6SlxR2, 
Therefore,  F  is  an  immersion  (6).  Clearly  F  is  injective. 
Since  IIF((e],  y,  z)ll2  =  e2y  +  z1  +  _y2  +  1.  we  have 


i Y  =  U-(?)  +£-K<<9)  (H2) 

where  Q  and  E=£0  are  constants.  Since  '.'*(*)  of  (17)  is 
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Fig  6  A  livallv  solvable  circuit  whose  inductor  constitutive  relation  is 
not  locally  current  controlled  (a)  The  circuit  diagram  (b)  Inductor 
constitutive  relation 


(a)  lh] 

Fig  ?  A  circuit  where  t,  is  not  a  (  ’  cinirdinale  la)  The  circuit 
diagram  lb)  Capacitor  constitutive  relation 


given  b\ 

SU) 

transversality  holds.  Since 


-M</  QY 


local  solvability  holds  and  the  dynamics  in  terms  of  q  is 
given  by 


^  =  -(q-Q)*-E.  (113) 


cause  ( ty  -  £ ),/J  +  Q  is  not  differentiable.  (It  is  not  even  a 
C'  diffeomorphism.)  On  the  other  hand  the  map: 

<£-('*•  vc.iK.i(..</) 

=  (-(</  (?)'t£. 

~{q~Q)'~  £•  (q  Q)'  E.q] 

is  a  C2  diffeomorphism  and  (113)  describes  the  dynamics. 
A  more  geometric  way  of  burking  at  the  situation  is  as 
follows.  The  configuration  space  2  is  diffeomorphic  to  the 
graph  of  g(  .  The  vector  field  Xiif  , ,  at  (Q.  E )  is  contained 
in  the  tangent  space  T(V  , ,2.  i.e..  XIV  , ,  is  in  parallel  with 
the  (/-axis.  Therefore,  if  we  look  at  A’,.,  , ,  from  the  </-axis. 
we  can  observe  the  direction  and  the  length  of  A'(,,  , ,.  On 
the  other  hand,  if  we  look  at  X(V  , ,  from  the  t(  -axis,  we 
can  detect  neither  the  direction  nor  the  length  of  A,,,  , 
Finally,  let  us  remark  that  even  though  (115)  does  not 
qualify  as  the  differential  equation  describing  the  dy¬ 
namics.  (115)  is  true  in  the  sense  that  for  the  flow  ,v(  t )  on 


<£<  (  v(t )) 
dt 


-  3(  t,  (  v(  t ))  E  Y  v(  i )). 


Our  geometric  approach  seems  to  be  the  right  tool  to 
explain  what  is  happening  in  this  example. 

We  will  next  give  a  perturbation  result  on  local  solvabil¬ 
ity.  Recall  that  \K  is  said  to  be  locally  hybrid  |3],  if  (9) 
holds  and 

det((  D/k  )A ), ,.k  t|( ,  for  all  (  ik )  S  A  K  i  L  K 

(11b) 


A  problem  arises,  however,  if  one  argues  as  follows:  Dif-  for  some  fixed  2n  K  *  n  K  matrix  A.  where  each  column  of  A 
ferentiating  (112)  with  respect  to  t  one  has  has  either  of  the  following  forms: 


=  l  D*c  Kj  =  ),  '<  =  (  *>Sch '« 

=  (/>*<•),«?*  =  -(*>*,  y),rt-.  (114) 

Since  =  '(i(  )  =  (t( ■-£)'  ’  +Q  and  since  (£>g(  )<<  = 

3 {q-Q)z.  one  has 

(*>*<•),  X,  o.  ,  =  -E\  \ 

This  and  ( 1 14)  imply 

~  =  -3(ty-£)2  'ty.  (115) 

Is  (115)  another  differential  equation  describing  the  dy¬ 
namics  of  the  same  circuit  as  ( 1 13)?  The  answer  is  no?  If  it 
were.  (113)  and  (115)  must  have  the  same  qualitative 
properties.  So  if  (115)  is  the  differential  equation  describ¬ 
ing  the  dynamics  of  the  circuit  of  Fig.  7,  then  t\  =  E  is  an 
equilibrium  of  the  dynamics.  But  the  corresponding  value 
q-%(  '(£)=()  is  not  an  equilibrium  point  of  (113).  Since 
the  existence  of  an  equilibrium  point  must  be  a  coordinate- 
free  property,  there  should  be  something  wrong  with  saying 
that  (115)  is  the  differential  equation  describing  the  dy¬ 
namics  of  the  circuit  of  Fig.  7.  The  point  here  is  that  ty  is 
not  a  C2  coordinate  for  2.  Although  the  map:  ty>-* 
<t>..  t(  •  if  <?)  =  (-  if  ty.  -  ty,  -  ty,(  ty  -  £)'  ’  +  Q) 
is  a  homeomorphism.  it  is  not  a  C2  diffeomorphism  be¬ 


ll).-  -.0.1.0.-  -.0.0. . 0 


(0. . .0.0.  -  -  .0. 1.0.  •  •  -  .0). 


Let 

*"7> 

and  suppose  that  f\  corresponds  to  i  H  (resp..  vH  ).  Then 
that  particular  resistor  is  said  to  be  locally  voltage  controlled 
(resp..  locally  current  controlled ). 

Remark:  Observe  that  in  (52)  and  (53).  local  controlled- 
ness  is  defined  for  \R.  whereas  in  the  above  definition, 
local  controlledness  is  defined  for  each  resistor  provided 
that  \R  is  locally  hybrid. 

Proposition  12:  Given  a  nonlinear  network  X ,  assume 
the  following. 

(i)  A*  is  locally  hybrid  and  A(  (resp.  A,)  is  locally 
charge  (resp.  flux)  controlled. 

(ii)  ADA>  0. 

Then,  by  adding  small  linear  capacitors  and  small  linear 
inductors  appropriately  we  can  obtain  a_  new  network  X 
such  that  (I)  ,\nAY  0.  (2)  A  ft  A,  (3)  X  is  locally  solva¬ 
ble. 

Proof:  Pick  a  proper  tree  V*  containing  a  maximum 
number  of  locally  voltage  controlled  resistors  and  a  mini- 
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mum  number  of  locally  current  controlled  resistors.  Let  i_ 
denote  its  associated  cotree.  Partition  (r.  i)  in  the  following 
manner: 


Hennrnts 

voltages 

cun  cuts 

localK  volume  controlled 

resistors  tn  r. 

*'*  t 

localK  current  controlled 

resistors  in  t: 

Jr 

inductors  in  c 

•i 

livdilv  voltage  controlled 

resistors  m  1 

«'u 

*i. 

locallv  current  controlled 

resistors  in '  > 

<-/, 

capacitors  in  ‘ ' 

r< 

it 

The  fundamental  loop  matrix  has  the  following  form: 


r, ,  r, 

«7  ‘V, 

‘7, 

*V 

1 

8// 

8/< 

:  «/. 

«/ 1 

8,, 

The  submatrix  B, ,  0  because  of  the  choice  of  the  tree. 

Now  insert  a  small  linear  capacitor  in  parallel  with  each 
locally  voltage  controlled  resistor  in  .'  and  insert  a  small 
linear  inductor  in  series  with  each  locally  current  con¬ 
trolled  resistor  in  t'.  Let  ~  /.,  U/,  U('U(',.  where  (.',  is 
the  branches  of  the  capacitors  added,  and  C  denotes  capa¬ 
citors.  It  is  clear  that  is  a  proper  tree  for  the  new 
network.  Statement  (1)  can  be  proved  in  a  similar  manner 
to  the  one  in  Proposition  4.  To  prove  (3)  observe  that  the 
fundamental  loop  matrix  for  -A  with  respect  to  V>  has  the 
following  form: 


1 


«•/, 


‘V, 


», 

8, 


«,< 

1 

8/(  8; , 


(m) 


where  /.,  represents  the  inductors  added.  Since  no  resistors 
are  added.  A  A  is  the  same  as  A*.  We  compute  the  matrix 
3CA(  t'A,  iA )  of  (100)  for  \A .  It  follows  from  ( I  P3)  that 


8, 


RR 


0.  Since  v„e  =(r,f.  tv, ).  p„,  =  ( r,, .  r,f ).  /, 


( /\e,  /, ,  )  and  .  i,e  )•  we  have 


l'/t  •  *«  )  .<,  )/«  (I-1 


20) 


Observe  that  the  matrix  of  (120)  depends  only  on  ( rA,  «„ ) 
and  that  if  ( vK,  iK )  6  irA(  £ ),  then  (rA,  i,)6A,.  Since  (120) 
is  obtained  simply  by  exchanging  columns  of  the  matrix  of 
(117).  it  follows  from  (116)  that  9C„(  t>.  i„ )  is  nonsingular. 
By  Corollary  7.  'A  is  locally  solvable.  This  proves  (3). 
Proposition  II  implies  AffiK  which  proves  (2).  U 

Example  II:  Consider  the  circuit  of  T.xample  2.  Since 
the  resistors  are  voltage  controlled,  insertion  of  and 
yields  local  solvability  of  the  circuit.  Therefore,  the  per¬ 


turbation  in  I  -.sample  2  was  already  good  enough  to  ensure 
not  only  Iransversality  but  also  local  solvability 

Example  12:  Consider  the  circuit  of  L. sample  3.  Insert¬ 
ing  ;i  small  linear  inductor  in  series,  one  can  make  the 
circuit  locally  solvable. 

Remark:  Ihe  number  of  reactive  elements  added  in 
Proposition  12  is  no  greater  than  the  number  of  reactive 
elements  added  in  Proposition  4.  Notice,  however,  that  in 
Proposition  12.  A A  is  required  to  be  locally  hvbrid.  whereas 
in  Proposition  4.  the  only  restriction  on  Aa  is  that  it  should 
be  an  n A -dimensional  C’  submanifold.  The  local  hybrid¬ 
ness  assumption  eannot  be  relaxed  as  the  following  exam¬ 
ple  shows. 

Example  Is:  Consider  the  circuit  of  big.  2(a)  where  the 
resistor  constitutive  relation  is  given  bv  the  unit  circle  S' 
(Lig.  8).  It  is  easy  to  check  A <Ti K.  In  fact  22  is  diffeomoi 
phic  to  S'.  Since  fH(  rA.  iK )  '  rK  I  and  since 

A  #(!/<.'«)  (  ^ihJr  hi  «.ikI  slp’L’  *)tat  del  ''■rI,r-ir)  ** 

at  points  .1  and  8.  Therefore,  the  circuit  is  not  locallv 
solvable.  Observe  that  A A  is  not  locally  hvbrid  since  there 
is  no  function  /A  satisfying  (^)  and  1 1 16)  for  a  fixed  .-I.  We 
claim  that  there  is  no  nut  of  making  the  circuit  locally 
solvable  bv  adding  linear  reactive  elements,  lo  show  this 
let  \  be  a  circuit  obtained  bv  adding  arbitrarv  number  of 
reactive  elements  to  the  original  circuit  ■\  .  Then  bv  ( 100). 
either  ’va(i:a.'k>  or  X  A<  f  A.  iH  ) 

(  /),  a  tH ),,  ^  ,  ,  depending  on  how  the  reactive  elements  are 
added.  In  any  case  there  are  points  where  del  'v  A.(  fA.iA  ) 

0.  'Therefore.  .\  cannot  be  locally  solvable. 

Note  that  the  perturbation  in  Proposition  12  is  a  network 
perturbation.  It  is  not  known  if  and  when  one  can  give 
element  perturbations  as  in  Proposition  s  in  such  a  muniiu 
that  ‘A  is  locallv  solvable.  One  can  xav.  however,  the  circuit 
of  bxaniple  13  cannot  be  made  locally  solvable  bv  element 
perturbations  To  see  this  let  a  perturbation  Aa  of  Aa.  .S'1 
be  described  by  /A(rA.iA)  0.  Since  S'  is  compact.  \A  is 
still  compact.  Therefore,  there  are  points  where 
(0.  >/«)„,,, .  0  and  «>.  Hence 

det  X  H{  rA.  i  A )  0  somewhere. 

If  A,  (resp.  A  /  )  is  not  locallv  charge  (resp.  flux)  con¬ 
trolled,  one  may  not  be  able  to  find  an  \  which  is  locallv 
solvable  as  the  following  example  shows. 

Example  14.  Consider  the  circuit  of  lag.  1(a)  where  the 
resistor  is  linear  and  A,  is  given  bv  lag  1(b).  Suppose  that 
linear  reactive  elements  are  added  in  such  a  wav  that  there 
is  still  a  proper  tree.  Let  (resp.  I. ,)  be  the  capacitance 
(resp.  inductance)  matrix  of  the  capacitors  (resp.  inductors) 
added.  Let  A  be  a  proper  tree  for  A  .  If  the  resistor  branch 
belongs  to  V>,  then  0C(  jr )  of  (S^)  is  given  bv 


*(*) 


:  «»/* 

f.  : 


Since  all  the  elements,  except  for  the  original  capacitor,  arc- 
linear.  and  since  \>  is  a  proper  tree,  one  can  show  that 
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w*(  22 )  and  for  (r.  1 ) C:  A'( r*.  /Jim1, 


rank 


<1 

Drfn  »./* 

On1'  />,«'' 


(  i:.m 


(r  ..rt  .  i,  )  serves  as  a  global  coordinate  system  for  — . 
Since  (  /)j»(  )f  l)  somewhere,  we  see  that  del  %( x  I  0 
there.  If  the  resistor  belongs  to  c\  then  a  similar  argument 
shows  that  vX  is  not  locally  solvable. 

We  will  next  discuss  relationship  between  local  solvabil¬ 
ity  of  X  and  transversality  of  the  (n<  t  n,  )-port  .V  derived 
from  'X  ,  under  certain  excitations.  Replace  capacitors  and 
inductors  of  'X  with  ports.  The  resulting  network  is  called 
the  («,  f  n,  )-port  V  derived  from  -X.  For  the  purpose  of 
convenience  we  will  keep  the  same  notation  for  .V  as  X . 
Drive  the  capacitor  ports  by  independent  voltage  sources 
r*  and  drive  the  inductor  ports  by  independent  current 
sources  1*.  Define 

•W;.!'*)-  {(r.i)eR:*i(«v#«)e.vt<  «;}• 

(121) 

This  set  represents  the  internal  constitutive  relations  of  V 
under  the  excitation  (r*.  1*).  Clearly  .V(  r*.  i* )  is  a  />- 
dimensional  submanifold.  Recall  22.  A  \  and  ir*  defined 
by  (72).  (74).  and  (76),  respectively. 

Proposition  13:  Given  a  nonlinear  network  \X .  assume 
that  A,  (resp.  A,  )  is  represented  by  q  gt  ( t;  )  (resp. 
<>  jf,(i(  ))  and  that  (Dgl  )r  (resp.  (Df>,  )i( )  is  symmetric 
and  positive  definite.  Assume  also  that  2i  is  an  (//,  t  n,  )- 
dimensional  C2  submanifold.  Then  ,X  is  locally  solvable  if 
and  only  if  for  the  <«,.  +  n,  )-port  ,V  derived  from  X.  the 
following  holds: 

At(f*./;)XAt.  for  all  (  r* ./"*)€**(—*)•  (122) 

Proof:  It  follows  front  the  hypothesis  and  Corollary  4 
that  \X  is  locally  solvable  if  and  only  if  ir*  is  a  local 
diffeomorphism  To  prove  sufficiency,  let  (r*.i*)€f,|-t) 
and  define 


By  an  argument  similar  to  that  of  the  proof  of  Proposition 
y.  one  can  show  that  ir*  is  a  local  diffeomorphism  if  and 
only  if  (125)  holds.  Finally,  since 

S'  A'nA'1  U  A'(  r*./*)i  A 1  (1261 

(I*. I?)'  n'l-'  I 

it  follows  from  (125)  that  for  each  22.  v'  is  a  local 

diffeomorphism.  Therefore.  X  is  Its.  alls  sols  able.  C  on¬ 
versely.  if  it*  is  a  local  diffeomorphism  at  each  (  r.  ip  22. 
then  (125)  and  (126)  imply  ( 122). 

V.  Fvi  nu  ai  Sirk  i  I’vssixm 

Fventual  strict  passivits  is  an  important  'juulitulisc 
property  of  electrical  networks,  because  it  guarantees  that 
all  trajectories  eventually  approach  a  fixed  compact  subset 
of  the  configuration  space  [12]  1 14).  Roughls  speaking,  the 
results  of  this  section  say  the  following:  Suppose  that  the 
resistors  are  eventually  strictly  passixe  and  that  esers 
capacitor  is  in  parallel  with  a  large  linear  resistor  and  esers 
inductor  is  in  series  w  ith  a  small  linear  resistor.  I  hen  all 
trajectories  approach  a  fixed  compact  subset  of  the  config¬ 
uration  space.  Since  the  above  assumption  is  satisfied  bs 
most  practical  networks,  the  results  guarantee  that  the 
voltage  and  current  waveforms  are  bounded  in  most  net¬ 
works  of  practical  interest. 

Consider  the  following  one-form  on  R  '  "■  : 

0  2  'v^t1  Sh.^i  »'-7> 

k  \  k  l 

and  suppose  that  capacitors  and  inductors  arc  reciprocal:4 

/?<  « 1 

de  £  <a-(  ♦  v  0  (i:s) 

k  I  k  I 


(i(  v.i)  - 


/r(  vk  •  'k  1 
«•*'(  e,  i) 


(125) 


where  w*':  R2*  -*R"'  1 is  the  projection  map 
ff*'(t\i)  -K ■.*/.) 


and  fH  is  defined  by  (9).  By  the  definition  of  A  1  ( r*.  /* ).  for 
each  ( r0.  »„  )G  .V(  v*.  1* ).  there  is  a  neighborhood  T'c  R :h 
of  ( q,. )  such  that 

A'(r;,i;)ni'  a  '(O) 

rank ( D(i )(rll  b.  for  all  ( r.  i ) G  A  '  ( r* .  ijf ) O  F. 

(124) 


Using  (122)  and  (121)  and  an  argument  similar  to  the 
proof  of  Proposition  I.  one  sees  that  for  each 


on  R:"“  x  A,  x  A, .  If.  in  addition.  21  is  simply  connected, 
then  [4)  there  is  a  unique  real-valued  function  E  on  21  such 
that  for  any  two  points  x  and  x„  of  2.’ 


where  T  is  any  smooth  curve  on  2i  connecting  x  and  jt„.  If 
we  fix  an  arbitrary  x„c2).  then  (129)  is  a  well-defined 
function  on  22  Clearly.  E  is  the  energy  stored  in  capacitors 
and  inductors  relative  to  the  point  x„  (E  2).  Let  If  j,:  R2"''  -» 
R  be  defined  by 

)*«( ~  2  (lit)) 

k  1 

4 Reciprocity  of  capacitors  and  induciors  is  related  wiih  the  existence  of 
energy  function,  whereas  reciprocity  of  a  network  1 10]  I  dec  0.  where  u  0 
defined  by  (00))  is  related  with  the  existence  of  mixed  potential  function 


4:: 
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fig  4  Diagram  defining  the  two  functions  ft  and  ft  „ 


Recall  vR  defined  by  (20)  and  let  tt'R  and  W  be  defined  by 
Fig.  9.  The  function  W  is  the  power  at  resistors.  It  follows 
from  Tellegen's  theorem  that 


</£'(*(*)) 

dt 


W'(x(f)). 


(131) 


Recall  that  a  network  '.X  is  said  to  be  eventually  strictly 
passive  [12]  [14|  if  there  is  a  compact  subset  such 

that 

fF(.r)>0.  forallxGS-fl.  (132) 

The  following  two  propositions  show  the  importance  of 
eventual  strict  passivity. 

Proposition  14:  \I2]-\I4)  Let  E  be  proper,  i.e.,  for  every 
n£R.  the  set  (x£-|  £(*)<<*}  is  compact,  and  let  ?X  be 
eventually  strictly  passive.  Then  the  set  defined  by 

P=  (x£2|  £(x)*£a,}  (133) 

a,  =  max£(x)  (134) 

at  St 

is  compact,  and  for  any  initial  state  x(0),  either  one  of  the 
following  happens: 

(i)  There  is  a  /,  >0  such  that  x(/)£f>  for  all  t>t{, 

(ii)  x(/)€c'>  for  all  f^O  but  lim,_xx(f)Efc\ 

The  set  t'-  contains  many  of  the  important  information 
concerning  the  dynamics.  In  particular,  the  following  holds: 

Proposition  15:  Under  the  same  setting  as  Proposition 
14.  we  have 

(i)  All  periodic  orbits  and  equilibria  are  in 

(ii)  In  particular,  equilibria  lie  in  the  set 

{xe2|ir(jr)  =  0}.  (135) 

Proof:  (i)  It  follows  from  (131)  and  (132)  that  for  any 
x(f)E2-ft.  the  energy  E(x(t ))  is  strictly  decreasing  with 
respect  to  t.  This  implies  that  for  x(0)E2-fi,  the  trajec¬ 
tory  x(t)  cannot  come  back  to  x(0).  Similarly,  x(t )  cannot 
remain  at  x(0). 

(ii)  Since  £(x(r))  is  either  strictly  increasing  or  strictly 
decreasing  outside  the  set  defined  by  (135).  the  equilibria 
must  be  located  in  ( 1 35).  □ 

The  set  f”  in  (133)  is  called  a  set  of  attraction  since  it 
attracts  all  trajectories. 

Eventual  strict  passivity  is  a  property  of  W  on  2.  while 
is  defined  on  A*.  These  two  functions  may  behave  very 
differently  depending  on  the  properties  of  t  and  vR.  (See 
Fig.  9).  The  properties  of  WR  are  much  easier  to  check  than 


Fig  10  A  network  which  is  not  eventually  strictly  passive  (a)  The 
circuit  diagram  (b)  Trajectories  on  the  linear  subspace  W  0 


those  of  W  because  W7,,  depends  only  on  \R  but  not  on  K 
so  that  one  does  not  have  to  worry  about  Kirchhoff  laws. 
We  need  the  following. 

Definition  4:  The  resistor  constitutive  relations  repre¬ 
sented  by  Ar  are  said  to  be  eventually  strictly  passive  if 
there  is  a  compact  subset  SlR  of  AR  such  that 

FF jt(  t-Ve  tR )>0,  for  all  ( vR .  iR ) £  A R  —  itR .  (13^, 

Eventual  strict  passivity  of  \R  is  a  physically  meaningful 
condition  because  it  simply  says  that  the  resistors  dissipate 
positive  power  eventually.  The  condition  is  satisfied  by  a 
broad  class  of  resistors.  A  natural  question,  then,  arises: 
Does  eventual  strict  passivity  of  A„  imply  existence  of  a 
compact  set  of  attraction?  Another  interesting  question 
related  to  this  one  was  raised  by  Smale  [8],  In  terms  of  our 
terminology,  the  problem  is  rephrased  as  follows:  Suppose 
that  there  is  a  number  fi>0  satisfying 

lVR(cR.iR)>0  2  K+dt)  (137) 

k  i 

for  all  (eR.  iR)  with  ll(t'R.  iR)ll  sufficiently  large.  Then,  does 
the  network  have  a  compact  set  of  attraction?  The  answer 
to  both  of  the  two  questions  is  no  as  the  following  example 
shows. 

Example  15:  Consider  the  circuit  of  Fig.  10(a).  where 
all  elements  are  linear  and  element  values  are  all  equal  to 
one.  Since  the  resistor  is  linear  and  1 -S2,  WR  is  positive 
everywhere  except  for  the  origin.  Hence  \R  is  eventually 
strictly  passive.  Observe  that 

^/*(  vr  -  'r  )  =  vr'r  =vR—'R^fl(vR  +  'r  ) 


for  0</3®S  I /2.  Therefore.  (137)  is  satisfied.  We  claim  that 
this  circuit  does  not  have  a  compact  set  of  attraction.  To 
this  end  let  us  write  the  dynamics  in  terms  of 
<‘V, .'i.): 


dlh 

dt 


d't.; 

dt 


«V,  *(■'#,  +  '/,)• 


(138) 


Drawing  trajectories,  (Fig.  10(b)).  one  can  show  that  there 
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A  R 

l  ig  1 1  Diagram  defining  the  t\u»  functions  %  and  i  K 

is  a  linear  subspace  on  which  all  the  trajectories  are  con¬ 
centric  circles.  More  specifically,  the  element  values  satisiv 
the  condition  for  this  bridge  circuit  to  be  balanced  at  the 
angular  frequency  one.  Therefore,  for  any  uOR.  the  fol¬ 
lowing  is  a  solution  to  ( 138): 

r ,,(/)  «',.(/)  asm  t 

i,  ( / )  if  j  t )  ucos/. 

Since  </£R  is  arbitrary,  the  solution  can  have  an  arbitrarily 
large  magnitude.  Therefore,  there  is  no  compact  set  of 
attraction.  In  terms  of  the  above  coordinate  system,  we 
have 

fT(  t<  .  i,  )  -  R(i,  .~i,  )' 

and  hence  it  does  not  satisfy  (132).  Notice  that  any  trajec¬ 
tory  starting  outside  the  linear  subspace  If  =  0.  approaches 
the  origin. 

Since  (13b)  is  satisfied  by  most  resistors  of  practical 
interest,  it  is  natural  for  us  to  seek  conditions  under  which 
(136)  implies  (132).  The  following  is  a  generalization  of  a 
recent  result  by  Chua  and  Green  [12]  for  a  general  mani¬ 
fold.  We  assume  that  \R  is  closed  for  technical  reasons. 
This  is  not  a  restrictive  condition,  however. 

Latum  2:  Let  21  be  an  (n,  +  n,  ^dimensional  (  ’  sub¬ 
manifold  and  assume  the  following: 

(i)  \R  is  closed  and  eventually  strictly  passive. 

(ii)  li(?.£>ll  -*  x  implies  !l(  r( .  i,  )ll  -t  x  on  At  v.\(. 
Then  -X  iy  eventually  strictly  passive  if  the  following 

fundamental  topological  hypothesis  is  satisfied: 

There  are  no  loops  and  no  cut  sets  consisting  only  of 
capacitors  and  inductors,  or  equivalently 

( 1 )  there  is  a  tree  (  R )  consisting  only  of  resistors. 

(2)  there  is  a  tree  ' (CL)  containing  all  capacitors  and 
inductors. 

Proof:  Recall  the  map  itR  defined  by  (22).  Suppose 
that  A„  is  eventually  strictly  passive  and  let  UR  be  as  in 
(136).  If  ttH  is  proper,  then  the  preimage  ttH  '(f2«)  is 
compact  because  the  preimage  of  a  compact  set  under  a 
proper  map  is  compact.  It  is,  then,  clear  that  the  inequality 
in  (132)  holds  with  respect  to  u,,  '($)„).  So  we  show  that 
the  fundamental  topological  hypothesis  implies  that  w#  is 
proper.  To  this  end  let 

•  (139) 

be  the  inclusion  map  and  consider  the  map  x  defined  bv 
Fig.  1 1.  Since  A  p  is  assumed  to  be  closed.  A  is  also  closed. 


4.’ I 

Therefore,  21  AO  A'  is  a  closed  submanifold  of  A.  Conse¬ 
quently.  for  any  compact  subset  -t  of  A.  the  preimage 

'(  A )  is  compact.  This  shows  that  iK  is  proper.  Therefore, 
we  need  only  show  that  \  is  proper.  Since  x  is  obviously 
continuous,  we  need  only  show  that  the  preimage  of  a 
bounded  subset  of  R:"*  is  bounded.  Suppose  that  the 
fundamental  topological  hypothesis  holds  and  let  (resp.. 
iKc )  be  the  tree  branch  voltages  (resp..  link  currents)  for 
'(  R)  (resp..  links  associated  with  ''((  /.)).  It  follows  from 
(15)  that  for  (  r.  i.if.  <j> )  L  A . 

«'  .  i  B'iR(  (140) 

where  Q  and  B  are  the  fundamental  cut  set  matrix  and  the 
fundamental  livop  matrix  associated  with  M  R  )  and  M  CL). 
respectively,  liquation  ( 140)  and  assumption  (ii)  imply  that 

!i.rll--x.  x  0  A  ->  !l(  *'«,  •  »*Kt )'!  — *  x  .  (141) 

Since  (  .  iR  )  is  a  subvector  of  (  rR.  iR  ).  we  have 

IIjtII— »oc.  xt  A ->  i(  cR.  iR  } 1  —  x .  (142) 

I  bis  shows  that  the  preimage  of  a  bounded  subset  under  x 
is  bounded.  Since  the  pioperties  of  x  do  not  depend  on  a 
particular  choice  of  a  tree,  x  is  proper. 

Remark:  Observe  that  in  the  above  proof  we  took  full 
advantage  of  the  coordinate  free  property,  since  in  (140) 
(142)  we  are  using  two  different  trees  simultaneously. 

Now.  experiences  tell  us  that  most  networks  of  practical 
interest  have  a  compact  set  of  attraction.  We  next  justify 
this  observation  formally  by  carrying  out  a  si  ght  network 
perturbation.  The  perturbation  we  make  is  simply  a  for¬ 
malization  of  the  following  hypothesis:  “livery  capacitor  is 
in  parallel  with  a  large  linear  resistor  and  every  inductor  is 
in  series  with  a  small  linear  resistor.”  Before  stating  the 
results,  we  need  the  following: 

Definition  5:  A  nonlinear  network  v\  is  said  to  be 
strongly  locally  solvable  if 

detOClx)^  0.  for  all  .rt  A  (143) 

where  DC(.tr)  is  defined  by  (K9)  and  A  is  defined  bv  (6). 

Remarks:  I)  If  A,  (resp.  A,)  is  locally  charge  (resp 
flux)  controlled,  then  "Ai  is  strongly  locally  solvable  if  and 
only  if 

detOCw(  vR,iR)=jtt(),  for  all  (  vR.  iR  )(_  \R  (144) 

where  0CW(  vR.  iR )  is  defined  by  ( 1(H)). 

2)  Condition  (143)  is  stronger  than  (k8)  since  for  strong 
local  solvability  the  determinant  should  be  nonzero  on  A 
and  since  2!cA.  This  condition  is  satisfied  bv  many  net¬ 
works.  however.  For  example,  the  circuit  of  Fig  2  with 
capacitors  added,  satisfies  this  condition  because 

„/k  |  ■ 

The  perturbed  network  X  of  Proposition  12  is  strongly 
locally  solvable  because  the  matrix  of  (120)  is  nonsingular 
for  all  ( t>.  iR )  £  A 

Proposition  16:  Given  a  nonlinear  network  A .  assume 
the  following: 

(i)  'X  is  strongly  locally  solvable. 
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(ii)  \R  is  closed  and  eventually  strictly  passive. 

(iii)  ll(fl.$)ll-»ao  implies  ||<t>(., »,  )ll -»  x  on  A,  x.\(. 
Insert  a  large  linear  resistor  gk ,  k=  1.-  •  -  ,«( ,  in  parallel 
with  every  capacitor  and  insert  a  small  linear  resistor  rk, 
k=  1.-  ■■.«,,  in  series  with  every  inductor.  Then 

(1)  A  O  K  0  and  2  =  A  D  K  is  an  («<+/;,)- 
dimensional  C*  submanifold, 

(2)  X  is  locally  solvable. 

(3)  X  is  eventually  strictly  passive.  Consequently  '.X  has 
a  compact  set  of  attraction. 

Proof:  (1)  In  terms  of  a  proper  tree  ‘5,  the  original 
network  '.X  is  described  by 

vRt  + Brrvr,  +  BrcV t  ~®  (145) 

V/  +  Bk RvR ^  q,  ;  0  ( 146) 

*r,  ~ —  B/r'i  =  ®  (147) 

«V (»4«) 

(t>.  i  ,q .  4>)  €.  A .  (149) 

Let  “T'-=  V'Ur,T  where  r.T  represents  branches  of  rk  s.  It  is 
clear  that  'X  is  a  proper  tree  for  .X .  Decompose  v  and  i  as 

c  =  (c^.cf,.c/:c/i,.cr.c(  ) 

H'/vV  ) 

where  g  and  /■  denote  the  variables  associated  with  gk's  and 
rk  s.  Then  X  is  described  by 


course.  ( 145)  ( 149)  are  written  as 
Br  If  v 

T  =0  (166) 
•  q  lin 

(v,  \.  (167) 

Comparing  ( 145 )- ( 149)  with  (150).  (153).  (I56).(I60) 
(167),  we  see  that  the  differences  between  X  and  X  are  in 
the  last  two  terms  of  ( 160).  ( 161 ).  and  (162)  ( 165).  There¬ 
fore.  'X  is  described  bv 


<?'  Jl'J  [G  J| 

(  l*’  Vr- ig'ir)~  H  ' 

*  1 


VR e  Vl  VKl  *< 


vRc  +  Brr vRf  +  Brcv(  -0 

(150) 

t',  +  B,  rvr^  +  BI  CV( -  +  cr  =  0 

(151) 

cs  +  c(  •  =  0 

(152) 

_  _  B[r*1  =  ® 

(153) 

h  Brc  *Re  —  Bi  k'i.  -**=® 

(154) 

»,-*/  ~  ® 

(155) 

(v.i,q.$)eA 

(156) 

'r=* 

(157) 

v,=rir 

(158) 

K=  diag(g,. •  •  ).  r=diag(r,.  -.rB)).  (159) 

Eliminating  rr  vr.  ig.  and  ir.  we  see  that  X  is  described  by 
(150).  (153),  (156).  and 

v,  +  B,rvr^  +  Blcv(  +  n,  =0  (160) 

>(~^RC*R,~^I.R*l.JtK  l,;<  =  ®  (161) 


Now  let  jc„  £  A  D  K  0  and  let  U  be  a  bounded  neighbor¬ 
hood  of  x{)  in  R 4 4"'.  Since  the  set  (AOA')ni'  is  a 
bounded  submanifold,  small  perturbations  of  K O  i  do  not 
destroy  transversality  of  A  n  U  and  A'OL'  and  hence  the\ 
do  not  destroy  nonemptiness  of  (  An  i '  )D(  AD  i  ).  There¬ 
fore.  if  gk  is  large  enough  and  if  rk  is  small  enough,  then 
H/XI  and  I1C II  in  (168)  are  small  enough  to  guarantee 
nonemptiness  of  the  intersection  of  (168)  and  (169).  Since 
(170)  does  not  destroy  this  nonemptiness.  we  have  A  n  A  =/- 
0.  We  show  AffiA  later. 

(2)  Since  (Dv  (resp.  (£>,,/)  )„•,.*,)  's  a  ^arc- 

matrix.  (89)  implies  that  if  X  is  strongly  locally  solvable, 
then  (0r,  A),,.,.  „,<resp.  (Di{  /,  )„;  *,)  is  nonsingular  on  A, 
(resp.  A, ),  i.e..  A(  (resp.  A, )  is  locally  charge  (resp.  flux) 
controlled.  Therefore,  we  can  use  Corollary  7  to  check  local 
solvability  of  X .  Observe  that 

/r(  vr  •  'R ) 

/*( •  >r  )  >„K  (173) 

t;  r  I, 

and  that  the  fundamental  loop  matrix  B  for  X  is  given  bv 


Let  us  rewrite  ( 145)-<  149)  and  ( 160)  ( 165)  more  concisely. 
Let  B  and  Q  be  the  fundamental  loop  matrix  and  the 
fundamental  cut  set  matrix  for  \X.  respectively.  Then,  of 


*>Re  vr 

«V  ""A 

V(  VR 
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It  follows  from  this  that 


7  V 

Y  R 

V 

\t2  / 

Bl  R-i^l  R  1  ]  Vt.  ■ 

Let  t’*e  =  (i V  V-  vR,  =(r«,.t-r).  iRe=0Rf,.ig).  and  i*, 
Then 

Di  „Jr  ~  (  ^i'kJr  )i*RR 


D*Jr 


D>hJr 


®r,,/«  (  D,  kJr  )Urr 


DiKefR  +  ( 

K,/. 


1>.»Jr+(D,-Jr)b£r 


It  follows  from  (176),  (177).  and  (100)  that 


Fig.  12.  Perturbation  of  network  of  Fig  I  (Ha). 

It  follows  from  condition  (ii)  that  there  is  a  compact  set 
fl^cAj  such  that  (136)  holds.  For  any  «''»().  let 


0ii/  .  . t 

Then  the  set  Q#  —  12 R  XS2xr  has  the  property  that 


WR  >0  on  A„  -ClR 

because  g  and  r  are  diagonal  matrices  with  positive  ele¬ 
ments.  Finally,  to  show  that  the  fundamental  topological 
hypothesis  is  satisfied,  let  Rg.  L.  r.-,.  and  g,  represent  the 
branches  of  the  resistors  in  ‘5.  inductors  in  L .  the  added 
resistors  rk 's  and  the  added  resistors  gk 's.  Then  (R)  =  R 
Ur.;  UgL  is  a  tree  for  which  consists  only  of  resistors. 
Also  X(CI.)  =  '5 U  L  is  a  tree  for  vK  which  contains  all 
capacitors  and  inductors.  It  follows  from  Lemma  2  and 
condition  (iii)  that  'X  is  eventually  strictly  passive  □ 

Example  16:  Consider  the  circuit  of  Example  15.  Since 
the  circuit  is  linear,  all  the  conditions  of  Proposition  1 6  are 
satisfied.  The  perturbed  circuit  is  shown  in  Fig.  12.  It 
follows  from  Proposition  16  that  this  perturbed  circuit  has  a 


^7? »  i  R  ) 

/«  -  (  dx  Jr  )brk  ■  [  d,jr  +  ( D,"  fR  )«’R  ■ 

1 

(17S) 

i  : 

It  is  clear  that  DC*  depends  only  on  ( vR.  iR )  and  that 

|det  9C„(  t>.  iK  )|  =  |det9C„(  vR.iR  )|.  (179) 

Now  if  ( i'R,  /„ )  E  A  then  ( r*.  iR  )£  A  K.  because 

*r  =  {(4-4)lbv4)e-V',  K  V rir). 

(ISO) 

By  the  strong  local  solvability  assumption,  we  have 
|det  0C#(  r* .  i„  )|  >0.  for  all  ( r „ .  iR )  £  A  „ . 

This  and  ( 179)  imply 

(del  0C*(  r*  ./„)!  >0.  for  all  ( r„ .  iR )  e  «„(  2 ). 

It  follows  from  Corollary  7  that  'T(  is  locally  solvable.  By 
Proposition  II.  we  have  A X .  Therefore,  2:  is  an  (n,  +n;  ) 
-dimensional  C2  submanifold. 

(3)  The  resistor  constitutive  relations  Afi  for  X  is  de¬ 
scribed  by  (ISO)  where  r*  (t>.  t'K.  tv).  iR  ( »«■  »,  )■ 

Therefore,  the  function  WR  corresponding  to  WR  is  given 

by 

iR  .iR)  H>(  vR  .iK)  +  v2  g  4  +  ijri, . 


compact  set  of  attraction.  In  fact  the  linear  subspace  If  "  0 
in  Fig.  10(b)  degenerates  into  the  origin  and  any  closed 
ball  centered  at  the  origin  serves  as  a  compact  set  of 
attraction. 

Remark:  As  we  have  seen  %„({■„.  iR)  %R(vR.iR)  for 
(e R.iR  )£«„(£).  But  (r„.  iR  )€«#(  —  )  docs  not  necessarily 
imply  (r s.i()6»#(i)  even  though  ( vR.  iM  )6  \H.  Recal¬ 
ling  Corollary  7  and  Definition  5.  one  sees  why  we  needed 
the  strong  local  solvability  hypothesis. 

We  will  next  replace  strong  local  solvability  with  another 
condition. 

Proposition  '7:  Replace  the  “strong  local  solvability” 
hypothesis  in  •'  roposinon  16  with  the  following  hypothesis: 

(i)’  v  is  a  global  diffeomorphism. 

Then,  under  the  same  perturbation  as  in  Proposition  16.  the 
same  conclusion  holds. 

Proof:  The  preceding  proof  for  Proposition  16  remains 
applicable  except  for  the  fact  that  £  is  an  ( n,  t  n,  )- 
dimensional  C2  submanifold  and  that  X  is  locally  solvable. 
In  order  to  prove  this  recall  (145)  (149).  Hypothesis  ( i )' 
implies  that  { v.  i.q.  <t>)  is  expressible  as  a  C2  function  of 
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Fig  13.  A  nonlinear  network  which  becomes  eventually  strictly  passive 
after  perturbation.  (a(  The  circuit  diagram  (b)  Resistor  constitutive 
relation. 


(«.<»); 

(v.i.q.  4t)  =  v~'  (q.<f>). 

Recall  ( 1 50>— ( 1 58)  and  set  v',  =  »,  +c,,  i'c  =  Then 

(150),  (151),  (153),  (154),  and  (156)  are  exactly  the  same  as 
( 145)— ( 149).  Therefore, 

K.%,«4-»*,io 4>)  =  w  '(q,4>).  (181) 

It  follows  from  (152),  (155),  (157),  (158)  and  (181)  that 
( r,,  vv  ij, )  is  also  a  C 2  function  of  (q.<t>).  Therefore,  all 
the  variables  of  X  are  expressible  as  a  C2  function  of 
(?.<>); 

(v.i.q.<t>)=v  ‘(q.<t>). 

If  follows  from  the  way  w  " 1  was  determined  that  it  1  is  a 
global  diffeomorphism  and  hence  so  is  v.  Therefore,  2  is 
an  ( nc  +  n ,  )-dimensional  C 2  submanifold.  Since  it  is  a 
global  diffeomorphism,  it  is_  a  local  diffeomorphism.  It 
follows  from  Theorem  I  that  \X  is  locally  solvable.  □ 

Example  17:  Consider  the  network  of  Fig.  13(a)  where 
the  resistor  is  described  by  Fig.  13(b).  The  resistor  is 
eventually  strictly  passive.  It  is  easy  to  show  that  w  is  a 
global  diffeomorphism.  Therefore  we  can  make  the  same 
perturbation  as  in  Example  16  so  that  the  network  will 
have  a  compact  set  of  attraction. 

We  will  next  relax  the  “strong  local  solvability”  hypothe¬ 
sis  and  the  global  diffeomorphism  assumption,  while  im¬ 
posing  a  stronger  condition  on  A,  to  derive  a  different 
perturbation  result.  Recall  that  AM  is  said  to  be  globally 
hybrid  [3]  if 

A*  =  {(»*.»*  ))/=*(  jt)} 

where  >’  =  (>v  ■  •■>'„„)•  *  =  (*,. •  •  „vn?)  and  if  t*  is  the 
current  (resp.,  voltage)  of  the  Ath  resistor  then  xA  is  the 
voltage  (resp.,  current)  of  the  Ath  resistor.  If  vk  is  the 
current  (resp.,  voltage),  then  that  particular  resistor  is 
called  voltage  controlled  (resp.,  current  controlled).  The 
following  result  says  that  most  practical  networks  can  be 
perturbed  in  such  a  manner  that  the  resulting  network  are 
locally  solvable  and  have  compact  sets  of  attraction. 

Theorem  2:  Given  a  nonlinear  network  :X  assume  the 
following: 

(i)  A„  is  closed,  globally  hybrid  and  eventually  strictly 
passive, 

(ii)  Ar  (resp.  A, )  is  locally  charge  (resp.  flux)  controlled 
and  ll(f. ♦)(!-* oo  implies  i,  )fl -» oo  on  A,  xA(, 

(iii)  An  0 . 

Perturb  X  in  the  following  manner: 

(a)  Let  5  be  a  proper  tree  containing  a  maximum  num¬ 
ber  of  voltage  controlled  resistors  and  a  minimum  number 


of  current  controlled  resistors  and  let  t  be  its  associated 
cotree.  Insert  a  small  linear  capacitor  in  parallel  with  each 
voltage  controlled  resistor  in  ‘,'i  and  insert  a  small  linear 
inductor  in  series  with  each  current  controlled  resistor  in  U. 
Call  the  resulting  network  'X . 

(b)  Insert  a  large  linear  resistor  gk  in  parallel  with  each 
capacitor  of  \X  and  insert  a  small  linear  resistor  rk  in  senes 
with  each  inductor  of  'X .  Call  the  resulting  network  'X . 
Then  the  following  hold:_ 

(1)  AnK^=0  and  2=  AHA  is  an  (n(  +n,  +  A  )- 
dimensional  C*  submanifold  where  A  is  the  number  of 
reactive  elements  added,  A  and  K  are  thej^esistor  constitu¬ 
tive  relations  and  the  Kirchhoff  space  of  -.X .  respectively, 

(2)  'X  is  locally  solvable. 

(3)  "X  is  eventually  strictly  passive.  Consequently  X  has 
a  compact  set  of  attraction. 

Proof:  (1)  It  is  clear  that  one  can  prove  AnA'^c  0  in 
a  similar  manner  to  the  proof  of  Proposition  16.  We  will 
prove  AffiA  later. 

(2).  (3)  We  first  claim  that  X  is  strongly  locally  solvable. 
To  this  end  partition  ( i\  i )  of  \X  as  in  the  proof  of 
Proposition  lb.  Since  A  R  is  assumed  to  be  globally  hybrid, 
it  can  be  represented  as  follows: 

'Y,  i/e)  =  ° 

)=0 

%  "0 

t7e”//f(t't,  l'lf  ® 

where  V  and  /  denote  voltage  controlled  and  current 
controlled  resistors,  respectively.  We  write  these  equations 
as 

/r(  vr  .  i  R)  ~  0. 

It  follows  from  (120)  that  for  X  we  have 


*Y  •  >R  )  -  D{v,,.r,t)fn  ;  4,', 


I <> • in  1 


Vl9  c/f  't  >v. 


y  1  ■  i 

Therefore.  3C„(  i'R.  iR)  is  a  constant  nonsingular  matrix  and 
"X  is  strongly  locally  solvable.  Clearly.  \R  =  Aj,  because  no 
resistors  are  added  in  (a).  This  implies  that  \R  is  eventu¬ 
ally  strictly  passive.  Since  ‘X  satisfies  the  hypothesis  of 
Proposition  lb.  by  taking  procedure  (b).  we  obtain  X  which 
is  locally  solvable,  A  A  and  eventually  strictly  passive.  □ 
Example  Id:  Consider  the  network  of  Fig.  14(a).  where 
and  R2  are  as  in  Fig.  2(b).  Other  elements  are  linear.  Bv 
a  similar  reasoning  to  that  of  Example  2.  one  can  show  that 
A  ftA.  Pick  the  proper  tree  X-  {C,.C2.C,,  R:).  Then 
applying  procedure  (a)  of  Theorem  2,  we  obtain  ‘X  which  is 
strongly  locally  solvable  (Fig.  14(b)).  The  network  X  of 
Fig.  14(b)  does  not  satisfy  the  fundamental  topological 
hypothesis,  however,  because  there  is  a  capacitor-only  cut 
set.  Insert  large  linear  resistors.  g2.  g,.  and  g4  according 
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Fig.  14  A  nonlinear  network  which  becomes  locally  solvable  and  even¬ 
tually  strictly  passive  after  perturbations  (a)  Original  network  X  (b) 
Perturbed  network  X  <c)  Perturbed  network  X 


to  procedure  (b)  and  obtain  vX  (Fig.  14(c)).  Theorem  2  says 
that  X  has  a  compact  set  of  attraction. 

Remark:  The  elements  added  in  Theorem  2  can  be 
thought  of  as  parasitic  elements  of  .  Therefore  Theorem 
2  formally  justifies  the  fact  that  in  most  networks  of 
practical  interest,  voltage  and  current  waveforms  eventu¬ 
ally  approach  a  fixed  compact  set. 

Appknuix 

Proof  of  Proposition  3:  We  will  first  prove  the  follow¬ 
ing: 

Lemma  A.  Suppose  that  A,  (resp.  A,  )  is  locally  voltage 
(resp.  current)  controlled.  Then  AffiX  if  and  only  if  A„ffi 
v'R(  K  ).  where  siR  is  defined  by  (20). 

Proof:  If  A(  (resp.  A, )  is  locally  voltage  (resp.  cur¬ 
rent)  controlled,  then  c(  (resp.  i, )  serves  as  a  local  coordi¬ 
nate  system  for  A(  (resp.  A,  ).  This  implies  that  (Dv(  )p 
(resp.  ( Di, )  )  in  (42)  is  nonsingular.  By  elementary  opera¬ 
tions,  one  can  show  that  (41)  holds  if  and  only  if 


Dl«e 

Brr 

Br( 

1 

MRt 

1 

p* 

o  r  of 

drh  °ir 

rank 


Next,  observe  that 


2/i , 


(A.l) 


*>*,  «V  'r{ 

~Brr  ~BR( 


1 

Bl,  B[r 

(A. 2) 


T{vK,iK  |Ar  -  Im 


Dv*i 

D>He 

Mr, 


(A. 3) 


Equations  (A. 2)  and  (A. 3)  imply  that  A R  ft  vR(  K )  if  and 
only  if  (A.l)  holds.  LJ 

Knowing  that  Affix  is  equivalent  to  A ffi  tt‘r{  X  ).  one 
sees  that  Proposition  3  can  be  proved  in  a  similar  manner 
to  that  of  Theorem  3  in  [3]  which  is  the  same  as  the  proof 
of  (ii-a)  of  Theorem  2  of  [3|.  Proof  of  (ii-a)  of  Theorem  2 
uses  Lemmas  /,  2,  and  4  of  (3J.  It  is  easy  to  show  that 
Lemma  !  is  true  for  C1  submanifolds.  Lemma  2  has 
nothing  to  do  with  differentiability.  Therefore  we  need  to 
only  show  that  Lemma  4  is  true  in  the  C2  category.  We 
state  this  in  the  following: 

Lemma  B:  Let  A  be  an  nXn  matrix  such  that  I \A  111 
is  sufficiently  small.  Then  there_are  neighborhoods  L\  and 
U2  of  the  origin  of  R"  with  (7,0(7,  and  there  is  a  C 2 
function  G:  R"  -»R"  such  that 

(i)  G-A  on  (7,, 

(ii)  G=id  off  (7,.  where  id  is  the  identity  map  of  R". 

(iii)  G  is  arbitrarily  close  to  id  in  the  strong  C 2  topology. 

Proof:  Let  %2(/<<;<(-))  be  a  sufficiently  small  neigh¬ 
borhood  of  ij  in  C2(R",R")  with  respect  to  the  strong  t 
topology.  Since  e(x)>0  for  all  x£R".  there  are  numbers 
r>0and  8>0  such  that  «(x)S*£  for  all  x  with  llx!l<8.  Let 
8„  satisfy  0<8o<8.  Then  there  is  a  C:  function  (bump 


function  [6))  p:  R"  — *R  such  that 

...  .  .  |l. 

*x) io.  UII/IIJ.S 

(A. 4) 

(ii)  There  is  a  k>0  such  that 

\\(Dp)x  II  <k,  \\(D2p)x\\<k 

(A. 5) 

for  all  xGR".  Now,  choose  A  close  enough  to  itl  so  that 


M-1IK 


t 

( J  4  8)(  I  +  2k~) 


(A. 6) 


and  define 

G(x)~  m(x)Ax  +  (I  -g(x|)x. 

We  will  show  that  G£‘}|2(i<j;  «( ■ )).  Since  g(x)  =  ()  for 
II  x  II  >8,  we  need  to  check  the  C2  size  of  G  id  only  for 
llxli<8.  Since  G(x)~x  =  p(xXAx  ~  x).  we  have,  using 
(A. 4)-  (A. 6),  that 


II  G(  x)  -  x  II  4  ||(  DG  )x  -  I II  +  ||(  DG  )i  II 
^h(x)IMx-x||  +  ||(Dm)xII  ll/lx-xll 
f  g(x)l(A  -  111  +  ||(Z)2ji)xll  IMx-xll  +  2II(  O/i  )x II 11/4  111 

<  11/4  -  111(11x11  +  k  llxll  +  I  +  k  llxll  +  2A  ) 

-  111(1  +8)(l  4  2A)<«. 

Take  (7,  -  {xeR',||lxll<6)  and  (A  (x€R"|  llxllcfi,,). 
Then  all  the  properties  are  satisfied.  □ 
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